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1.0  SUMMARY 


An  experimental  program  was  conducted  to  investigate  heat  transfer  characteris¬ 
tics  of  rotating  multipass  passages  for  configurations  and  flow  pareimeters 
typical  of  modem  gas  turbine  blades.  The  experiments  were  conducted  with  a 
smooth  wall,  large  scale  heat  transfer  model.  The  objective  was  to  obtain  the 
heat  transfer  data  base  required  to  develop  heat  transfer  correlations  and  to 
assess  computational  fluid  dynamic  techniques  for  rotating  coolant  passages. 

An  analysis  of  the  governing  equations  showed  that  four  parameters  influence 
the  heat  transfer  in  rotating  passages  (coolant  density  ratio,  rotation  number, 
Reynolds  number  and  radius  ratio).  These  four  parameters  were  varied  over 
ranges  which  exceed  the  ranges  of  current  open  literature  results,  but  which 
are  typical  of  current  and  advanced  gas  turbine  engine  operating  conditions. 
Rotation  affected  the  heat  transfer  coefficients  differently  for  different 
locations  in  the  coolant  passage.  For  example,  heat  transfer  at  some  locations 
increased  with  rotation,  but  decreased  and  then  increased  again  at  other 
locations.  Heat  transfer  coefficients  varied  by  as  much  as  a  factor  of  5 
between  the  leading  and  trailing  surfaces  for  the  same  test  condition  and 
streamwise  location.  Comparisons  with  previous  results  are  also  presented. 

This  work  was  supported  by  the  NASA/Lewis  Research  Center  under  the  Hot  Section 
Technology  (HOST)  initiative,  Contract  No.  NAS3-23691  to  Pratt  &  Whitney, 
Commercial  Engineering.  The  work  was  performed  under  the  direction  of  Mr.  F. 
Yeh,  NASA  Project  Manager,  and  Mr.  S.  Tanrikut,  Pratt  &  Whitney  Progr£un 
Manager.  Additional  experiments  were  conducted  under  United  Technologies 
Corporation  independent  research  sponsorship  in  order  to  enhance  the  benchmark 
data  base.  These  data  are  also  included  in  this  report. 
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2.0  INTRODUCTION 


Current  high  performance  gas  turbine  engines  exploit  internal  convection 
cooling  schemes  to  maintain  acceptable  airfoil  metal  temperatures.  This 
approach  requires  complex  coolant  passage  configurations  within  the  rotating 
blades  as  Illustrated  in  Figure  1.  Most  coolant  passage  designs  enhance  heat 
transfer  coefficients  above  smooth  channel  levels  by  utilizing  turbulence 
promoters.  Summaries  of  the  technical  problems  and  needs  of  the  current  gas 
turbine  blade  designer  are  presented  by  Suo  (1978)  and  Taylor  (1980). 

A  review  of  general  heat  transfer  literature  by  Bergles  and  Webb  (1970), 
contains  significant  references  for  turbine  airfoil  cooling  passage  design. 

Webb  et.al.  (1971)  has  reported  measurements  and  correlations  for  flow  in 
nonrotating  tubes  with  repeated  rib  roughness.  Burggraf  (1970)  and  Han  et.al. 
(1978)  conducted  experimental  heat  transfer  studies  with  rib  roughened 
geometries  typical  of  gas  turbine  engines.  Buoyancy  effects  in  vertical 
stationary  ducts  were  reported  by  Eckert  et.al.  (1953),  Metals  and  Eckert 
(1964)  and  Brundrett  and  Burroughs  (1967).  There  are  limited  amounts  of 
rotating  passage  heat  transfer  data,  with  the  bulk  of  this  work  done  with 
circular  tubes.  The  effects  of  rotation  on  secondary  flow  and  stability  have 
been  investigated  by  Moore  (1967),  Hart  (1971) , Wagner  and  Velkoff  (1972), 
Johnston  et.al.  (1972)  and  Rothe  and  Johnston  (1979).  Heat  transfer  in 
rotating,  smooth  wall  models  has  been  investigated  by  Mori  et.al.  (1971), 
Johnson  (1978),  Morris  and  Ayhan  (1979),  Lokai  and  Gunchenko  (1979),  Morris 
(1981),  Iskakov  and  Trushin  (1983)  and  more  recently,  Guidez  (1988).  Some  of 
the  results  contained  in  this  report  have  been  previously  reported  in  Wagner 
et.  al.  (1989,  1991)  and  Wagner  et.  al.  (1990). 

Large  Increases  and  decreases  in  local  heat  transfer  were  found  to  occur  by 
some  Investigators  under  certain  conditons  of  rotation  while  others  showed 
lesser  effects.  Analysis  of  these  results  does  not  produce  consistent  trends 
in  the  effects  of  rotation  on  heat  transfer.  The  disparity  in  the  results  is 
indicative  of  differences  in  the  measurement  techniques  and  models  used  in  the 
experiments  as  well  as  the  nonuniformity  of  the  test  conditions. 

Initial  work  to  investigate  rotating  boimd  shear  flows  ex2unined  unheated 
circular  and  rectangular  duct  flows.  For  radially  outward  flow,  rotation  was 
shown  to  generate  secondary  flows  perpendicular  to  the  mainstream  flow 
direction  in  ducts  orthogonal  to  the  axis  of  rotation.  Much  of  the  early  work 
was  conducted  for  laminar  flow  since  rig  limitations  in  an  unpressurized  duct 
allow  only  low  Reynolds  nvunber  flows  at  appropriate  rotation  or  Rossby  numbers. 
The  secondary  flow  patterns  of  Barua  (1954-1955),  shown  in  Figure  2,  illustrate 
the  effect  of  rotation  on  isothermal  laminar  flow  in  a  circular  duct.  Ito  and 
Nanbu  (1970)  conducted  one  of  the  first  heat  and  mass  transfer  experiments  in 
a  rotating  circular  tube.  These  results,  shown  in  Figures  3  and  4,  compare  the 
non-rotating  and  rotating  values  of  local  heat  transfer.  The  small  dashed  inner 
circle  in  Figure  4  represents  the  level  of  non-rotating  heat  transfer.  Notice 
the  strong  shift  in  heat  transfer  in  the  rotating  duct.  Plotted  on  a  polar 
plot,  most  of  the  tube  exhibits  significant  Increases  in  local  Nusselt  nvunber 
up  to  six  times  the  non-rotating  level.  Only  precisely  on  the  centerline 
of  the  leading  surface  is  a  decrease  in  heat  transfer  evident. 
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Figure  1.-  Cooling  Concept  of  a  Modern  Multipass  Turbine  Blade. 
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Figure  2.-  Typical  Secondary  Flow  Pattern  For  Isothermal  Laminar  Developed 
Flow  With  Orthogonal-Mode  Pattern  Rotation  (Barua,  1955). 
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Figure  3.-  Local  Mass  Transfer  Result  in  Laminar  Region  (Ito  and  Nanbu,  1970). 
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Figure  4.-  Local  Analogous  Nusselt  Nximber  in  Laminar  Regions  (rotation) 
(Ito  and  Nanbu,  1970). 


The  generation  of  secondary  flows  by  Coriolis  forces  produces  a  complicated 
three-dimensional  flow  that  greatly  alters  the  circiamferential  heat  transfer 
distribution  in  rotating  duct  flows.  The  secondary  flows  produced  in  these 
configurations  will  influence  the  distribution  of  heat  transfer  according  to 
the  movement  or  migration  of  fluid  from  one  surface  to  another  and  the  mixing 
of  the  near-wall  fluid  with  the  mainstream  or  core  flow.  Thus,  the  circumferen¬ 
tial  distribution  of  heat  transfer  will  be  vastly  different  in  each  of  these 
geometries,  as  it  will  be  in  the  wide  variety  of  cooling  channel  geometries 
found  in  advanced  aircraft  gas  turbine  blades. 

Most  of  the  published  work  on  rotating  duct  flows  has  been  on  experiments 
conducted  with  long  straight  ducts  with  flow  traveling  either  radially  outward 
or  radially  inward  with  respect  to  the  axis  of  rotation.  As  depicted  in  Figure 
1,  modem  gas  turbine  airfoils  employ  complicated  serpentine  shaped  coolant 
passages  that  utilize  coolant  flowing  both  radially  Inward  and  radially 
outward.  During  the  rotation  of  blades,  the  centrifugal  forces  of  rotation 
will  influence  the  heat  transfer  in  each  of  these  cases  differently.  Also,  the 
upstream  and  downstream  influence  of  the  turns  will  affect  the  heat  transfer 
in  these  passages  and  becomes  very  complicated  during  rotation.  A  further 
discussion  of  the  physical  effects  of  rotation  on  heat  transfer  including  the 
work  of  other  experimenters  is  presented  in  Sections  8.1  and  8.3. 
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Currently,  the  analysis  of  airfoil  internal  passage  heat  transfer  and  pressure 
loss  relies  mainly  on  correlations  derived  from  testing  models  in  a  static 
(nonrotating)  environment.  Executing  tests  with  rotation  is  difficult  and 
costly.  As  a  consequence,  there  is  limited  data  that  can  be  used  to  account 
for  the  effects  of  rotation  on  internal  heat  transfer  and  pressure  loss  in 
typical  turbine  blade  designs.  Some  data  are  available  for  smooth  tubes  over  a 
limited  range  of  revelant  parameters,  but  application  of  these  data  to  compli¬ 
cated  flow  passages  of  a  turbine  airfoil  would  not  be  appropriate.  Presently, 
adjustment  factors  are  applied  to  the  static  test  derived  correlations  to  bring 
them  into  nominal  correspondence  with  engine  experience.  This,  in  practice, 
accounts  for  rotation  effects. 

The  objective  of  this  phase  of  the  program  was  to  acquire  heat  transfer  and 
pressure  drop  data  and  to  develop  correlations  for  multipass  rotating, 
smooth-surface  coolant  passages  under  conditions  similar  to  those  expected  in 
the  first  stages  of  advanced  aircraft  gas  turbines.  Local  heat  transfer  was 
measured  along  the  smooth  wall  coolant  passage  and  around  its  periphery  for 
radial  outflow  and  inflow  conditions.  Local  heat  transfer  was  also  measured  in 
the  three  turns,  two  at  the  model  tip  and  one  at  the  model  root.  Incremental 
pressure  drop  was  measured  along  the  passages.  All  data  were  obtained  for  a 
range  of  rotational,  flow,  and  heat  fltuc  conditions  representative  of  engine 
operation. 

The  information  generated  with  the  smooth  model  was  evaluated  and  compared 
with  the  results  from  the  following  two  test  phases  utilizing  turbulators  in 
the  same  model  to  enhance  local  heat  transfer.  Upon  completion  of  the  entire 
program,  detailed  information  will  be  available  on  the  effects  of  rotation  on 
internal  flow  and  heat  transfer  in  rotating,  heated  passages  with  and  without 
turbulators.  A  complete  set  of  data  files  from  all  three  phases  of  this 
experimental  prograun  can  be  obtained  through  Mr’.  F.  Yeh,  NASA  Project  Manager 
at  the  NASA  Lewis  Research  Center. 
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3.0  EXPERIMENTAL  APPARATUS  AND  PROCEDURES 


3.1  Heat  Transfer  Model 

A  four  legged  heat  transfer  model,  used  for  the  NASA-sponsored  heat  transfer 
and  pressure  drop  experiments,  was  designed,  fabricated  and  Instrument a ted 
under  UTC  sponsorship.  The  model  consists  of  three  heated  straight  sections, 
one  unheated  straight  section,  and  three  turn  sections  as  shown  In  Figure  5. 
The  streeunwlse  location  of  each  test  surface  Is  Identified  by  a  letter  A  to  R. 
The  orientations  of  the  test  surface  at  each  streamwlse  location  are  denoted 
"leading"  and  "trailing"  for  the  surfaces  In  the  plane  of  Figure  5  and  "side 
walls"  for  the  surfaces  (crosshatched)  perpendicular  to  the  plane  of  Figure  5. 
Cross  sections  of  the  straight  sections  are  shown  In  Figure  6.  The  model 
was  designed  for  constant  temperature,  steady-state  heat  balance  measurements 
and  for  wall  static  pressure  measurements.  A  photograph  of  the  uninstrumented 
coolant  passage  heat  transfer  model  Is  shown  In  Figure  7. 

3.1.1  Model  Design 

Each  streamwlse  location  In  the  straight  heated  sections  has  the  cross- 
sectional  shape  and  features  shown  In  Figure  6.  All  four  copper  walls  were 
heated  on  the  side  opposite  the  test  surface  with  thin  film  electric 
resistance  heaters.  The  heaters  were  designed  to  produce  a  maxlmiun  heat  flux 
of  4.6  to  6.2  watts  per  square  centimeter  (30  to  40  watts  per  square  Inch). 

The  heaters  were  fastened  to  the  copper  test  surfaces  using  standard  strain 
gage  adhesives  useable  for  177®C  (350®F)  operation.  The  temperatures  of  the 
copper  test  surfaces  were  measured  with  two  chromel-alumel  thermocouples 
Inserted  Into  drilled  holes  of  each  test  surface  and  fastened  with  epoxy.  The 
copper  test  surfaces  were  separated  from  each  other  In  both  the  streamwlse 
direction  and  around  the  coolant  passage  cross  section  with  1.52mm  (0.060  In.) 
thick  sections  of  6-10  or  G-11  leunlnated  fiberglass  material.  Details  on  the 
heat  balance  and  the  calculation  of  the  effective  heat  transfer  area  for  each 
test  section  will  be  discussed  In  a  subsequent  section. 


The  turn  sections  had  three  sides  of  the  turn  cross  section  heated  and  the 
fourth  side  unheated.  The  unheated  surface  was  the  Inner-radius,  curved 
surface  shown  In  Figure  5.  One  heated  test  surface  (leading  and  trailing  wall 
surfaces)  covers  the  coolant  passage  In  the  plane  shown  In  Figure  5  for 
streamwsle  locations  E,  F,  J,  K,  P  and  R.  Two  heated  test  surfaces  cover  the 
outer  radius  curved  surfaces  (side  walls  surfaces). 
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TEST  SECTION  ELEMENT  IDENTIFICATION 

SIDE  WALL  TEST  SECTION  SURFACES.  1-32  ARE  IN 

PLANE  PERPENDICULAR  TO  VIEW  SHOWN 

TEST  SECTION  SURFACES  33-48  ARE  ON  "  +  Q"  LEADING  PLANE 
TEST  SECTION  SURFACES  (49)-(64)  ARE  ON  ”-1-  2  '  TRAILING  PLANE 
PRESSURE  MEASUREMENT  LOCATIONS^  -  fe] 


Figure  5.-  Cross  Sectional  View  of  Coolant  Passage  Heat  Transfer  Model 

Assembly.  Streamwise  Location  of  Test  Sections  Identified  by  A  to 
R.  All  Four  Test  Section  Surfaces  for  Streamwise  Locations  A 
Through  R  are  Heated. 
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Figure  7.-  Photograph  of  Uninstrumented  Coolant  Passage  Heat  Transfer  Model 
With  Leading  Edge  (+n)  Plane  Test  Sections  Removed. 

t 

3.1.2  Model  Instrumentation 

The  instrumentation  plan  and  the  test  surface  instrumentation  plan  for  the 
coolant  passage  heat  transfer  model  are  shown  in  Figures  8  and  9,  respectively. 
The  streamwise  locations  (A  through  R)  previously  shown  in  Figure  5  are  repeat¬ 
ed  in  Figure  8.  In  addition  each  heat  transfer  test  surface  is  identified  by  a 
nximber  (1  through  64).  Note  that  test  surfaces  2,  20,  34  and  50  are  located  at 
streamwise  location  B  (Figure  8).  Test  surfaces  34  and  50  are  the  leading  and 
trailing  surfaces  respectively  when  the  model  is  in  the  Q!=  0®  orientation 
(Figure  9).  Test  surfaces  2  and  20  are  at  sidewall  locations  for  a  =  o®. 

Pressure  measurement  locations  are  shown  at  16  places  on  Figure  8.  The  loca¬ 
tions  are  on  the  side  walls  as  shown  and  were  chosen  to  separate  the  turn 
pressure  losses  from  the  straight  section  losses.  The  pressure  tap  is  a 
0.51  mm  (0.020  in.)  diameter  hole  drilled  perpendicular  to  the  test  section 
surface  midway  across  the  passage. 
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TEST  SECTION  ELEMENT  IDENTIFICATION 

SIDE  WALL  TEST  SECTION  SURFACES.  1-32  ARE  IN 

PLANE  PERPENDICULAR  TO  VIEW  SHOWN 

TEST  SECTION  SURFACES  33-48  ARE  ON  "-i-Q''  LEADING  PLANE 
TEST  SECTION  SURFACES  (49)-(64)  ARE  ON  " -f  Q  '  TRAILING  PLANE 
PRESSURE  MEASUREMENT  LOCATIONsjT]  -  jfs] 


LOCATION  OF  THE 
END  OF  A  SQUARE 
CHANNEL  MOUNTED 
ON  THE 
INLET  NOZZLE 


HEATED  STRAIGHT 
TEST  SECTIONS 


EXIT  OF  THE 
INLET  NOZZLE 
SYSTEM 


EATED 
SECTIONS 


MEASURED  INLET 
BULK  TEMPERATURE 


MEASURED  OUTLET 
BULK  TEMPERATURE 


Figure  8.-  Instrumentation  Plan  for  Coolant  Passage  Heat  Transfer  Model. 


11 


a  =  0  DEG.  ORIENTATION 


FIRST  LEG 
FLOW  OUTWARD 


SECOND  LEG 
FLOW  INWARD 


THIRD  LEG 
FLOW  OUTWARD 


FOURTH  LEG 
FLOW  INWARD 


33-37 


38-42 


ROTATION 


019-22_rP^8J1_ 

23-26Lj^ _ 


43-47 

=P 


48 


49-53 


V 


54-58 


TEST 

SECTION 

SURFACE 


59-63 

AXIS  OF  ROTATION 


Figure  9.-  Test  Surface  Identification  Plan  for  Coolant  Passage  Heat  Transfer 
Model  (View  from  the  axis  of  rotation  looking  radially  outward  to 
the  model) . 
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The  test  section  surface  locations  for  the  a  =  0®  and  the  a  =  45®  orientations 
are  shown  in  Figure  9.  Note  that  for  0!=  0®,  the  centerlines  of  all  four 
straight  test  sections  pass  through  the  axis  of  rotation  and  the  straight 
sections  are  also  radial  sections.  For  Oi  =  45®,  the  centers  of  the  test 
sections  are  offset  from  radial  as  shown  in  Figure  9.  The  radial  positions  of 
several  locations  are  shown  on  Figure  5  for  reference. 

A  photograph  of  the  partially-assembled,  instrumentated,  coolant  passage  heat 
transfer  model  is  shown  in  Figure  10.  Each  of  the  64  test  surfaces  has  two 
chrorael-alumel  thermocouples  imbedded  and  has  a  thin  film  heater  attached.  The 
six  leads  from  each  test  surface  multiplied  by  sixty  four  test  surfaces  result 
in  384  leads  for  this  portion  of  the  model.  Additional  thermocouples  are 
positioned  at  the  air  inlet  and  exhaust  locations  and  on  the  steel  backing 
plates.  A  photograph  of  the  assembled  model  mounted  on  the  base  is  shown  in 
Figure  11.  A  photograph  of  the  model  mounted  in  the  Rotating  Heat  Transfer 
Facility  is  shown  in  Figure  12. 

3.2  Rotating  Heat  Transfer  Facility 

The  Rotating  Heat  Transfer  Facility  (RHTF)  (Figure  13)  consists  of  the 
containment  vessel  with  the  Integral  arm  assembly  and  motor  with  associated 
controller.  The  containment  vessel  is  1.83  m  (6.0  ft.)  in  diameter  and  was 
designed  to  withstand  a  destructive  failure  of  the  rotating  assembly.  The 
vessel  was  designed  for  operation  at  pressure  of  5  to  13  mm  of  Hg  absolute  to 
reduce  the  power  required  to  rotate  the  arm.  The  rotating  arm  assembly  is 
driven  by  a  IIKW  (15  Hp)  DC  motor  via  a  toothed  belt.  Shaft  rpm  is  controlled 
by  an  adjustable  feedback  electronic  controller.  Maximum  shaft  speed  is 
approximately  3,500  rpm  producing  body  forces  on  the  model  of  approximately 
14,000  g's  at  the  tip  of  the  model  and  approximately  10,000  g's  at  the  root. 
The  maximum  shaft  speed  for  the  present  program  was  1100  rpm.  A  safety 
shutdown  interlock  circuit  is  used  to  turn  off  the  drive  motor  and  model 
heater  power  supplies,  turn  on  a  magnetic  brake  and  open  the  containment 
vessel  vacuum  chamber  vent.  The  safety  shutdown  system  prevents  damage  to  the 
model  or  the  facility  in  the  event  of  a  leak  in  the  model  or  an  imbalance  in 
the  rotating  assembly. 

The  shaft  assembly  comprises  a  main  outer  shaft  with  two  shorter  inner  shafts. 
This  shaft  arrangement  was  designed  for  dual  fluid  paths  from  each  rotary 
union  mounted  on  the  ends  of  the  shaft  to  the  rotating  assembly.  Grooves 
located  on  the  exterior  surface  of  the  outer  shaft  allow  instrumentation  and 
power  leads  to  extend  from  the  rotating  arm  to  the  rotating  portion  of  the 
instrumentation  slipring.  Two  slipring  assemblies  (a  40  channel  unit  located 
on  the  upper  end  of  the  shaft  and  a  200  channel  unit  located  on  the  lower  end 
of  the  shaft)  are  used  to  transfer  heater  power  and  instrumentation  leads 
between  the  stationary  and  rotating  frames  of  reference. 
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Figure  10.-  Photograph  of  Instrumented  Coolant  Passage  Heat  Transfer  Model 
With  Trailing  Edge  (+0)  Plane  Test  Section  Removed  (Two 
thermocouples  and  thin  film  heaters  mounted  on  each  test  section) 


ELECTRICAL 
CONNECTORS 
(2  TYPICAL! 


PNEUMATIC 
CONNECTOR  (11 


BASE 


Figure  11,-  Photograph  of  Assembled  Model  Mounted  on  Base  With  Pressure  Shell 
Removed. 


Figure  12.-  Photograph  of  Model  Mounted  in  Rotating  Heat  Transfer  Facility 
(Rotating  heat  transfer  facility  with  cover  removed). 


Figure  13,-  Rotating  Heat  Transfer  Facility. 


3.2.1  Data  AcQulsition  System 


The  data  acquisition  system  contains  two  major  components;  the  computer  and 
the  data  acquisition  control  unit.  The  computer  consists  of  a  DEC  PDF  11/03 
processor  unit  with  128k  memory,  two  20cm  (8  in.)  floppy  disk  drives  and  a 
DECWRITER  III  terminal.  The  Hewlett  Packard  3497A  data  acquisition  system  can 
be  controlled  from  the  front  panel  or  through  the  interface  connected  to  the 
computer.  The  model  heater  voltages  are  set  manually  and  adjusted  until  the 
required  wall  surface  temperatures  are  obtained.  Upon  completion  of  the  acqui¬ 
sition  of  voltage  data,  results  are  calculated  and  printed  in  engineering 
units.  Flow  parameter  and  raw  data  are  stored  on  disk  for  future  reduction. 

3.2.2  Heater  Power  Source 


The  power  supply  system  provides  DC  power  for  the  thin  film  foil  resistance 
heaters  used  to  heat  the  model  test  section  elements.  There  are  72  individually 
controlled  power  supplies  which  are  rated  for  50  watts  of  power  with  a  meucimum 
current  draw  of  one  amp.  Individual  units  can  be  arranged  in  parallel  as  needed 
to  supply  additional  power.  Heater  supply  voltage  and  the  voltage  across  preci¬ 
sion  current  measurement  resistors  are  measured  by  the  data  acquisition  unit. 

3.2.3  Flow  Monitoring  System 

Model  coolant  air  is  supplied  by  the  UTRC  27  atm  (400  psig)  air  system  which 
is  regulated  to  approximately  10  atm  (150  psig)  at  the  RHTF.  The  air  flow  rate 
is  measured  with  variable  area  flow  meters.  The  model  coolant  return  air  flows 
through  an  additional  flow  meter  to  determine  a  mass  flow  balance  on  the 
system.  Model  pressure  is  controlled  by  back  pressuring  the  model  air  flow 
system  with  the  return  air  control  valve.  The  maximum  mass  flow  rate  available 
is  dependent  on  the  model  operating  pressure  and  the  total  pressure  loss  of 
the  system  including  the  heat  transfer  model.  For  typical  models  the  maximtun 
air  flow  rate  is  approximately  0.02  kg/sec.  (0.044  Ibm/sec) 

I 

3.3  Experimental  Procedures 

Heat  transfer  characteristics  was  determined  from  a  heat  balance  on  each 
heated  test  surface.  The  heat  added  to  the  coolant  by  convection  was  deter¬ 
mined  from  the  electrical  power  used  to  heat  each  test  surface  and  the  heat 
conducted  from  the  test  surface  to  the  support  structure.  Heater  supply  power 
for  individual  model  segments  was  determined  by  multiplying  the  calculated 
voltage  across  each  of  the  microfoil  heater  leads  and  the  current  determined 
from  the  voltage  measured  across  precision  0.1  ohm  resistors.  The  voltage 
across  the  microfoil  heater  leads  was  determined  by  accounting  for  the  voltage 
drop  across  the  heater  supply  leads.  The  net  heat  flux  (convected  heat  flux), 
assuming  negligible  radiation  energy  transfer,  was  determined  by  subtracting 
the  conducted  backloss  from  the  heater  power  input. 
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The  conduction  loss  parameter  for  each  heated  segment  was  determined  by  uni¬ 
formly  heating  the  model  segments  with  no  coolant  flow  to  a  steady  state 
condition  and  measuring  the  voltage  and  current  required  to  calculate  heater 
power  and  all  the  model  temperatures.  For  this  condition,  the  convected  heat 
flux  is  zero  and  the  total  heater  power  is  backloss.  The  conduction  backloss 
parameter  is  calculated  by  dividing  this  heater  power  by  the  temperature 
difference  of  the  heated  segments  and  the  support  frame.  The  bulk  temperature 
used  to  calculate  heat  transfer  coefficients  was  determined  with  a  thermo¬ 
dynamic  energy  balance  through  each  discrete  system  of  heated  segments. 
Rotating  heat  transfer  results  were  normalized  with  the  nonrotating  values. 

Electronic  noise  in  the  data  signals  of  the  RHTF  was  present  only  when  the 
shaft  was  rotating.  The  probable  sources  for  this  electronic  noise  were  (1) 
rotating  instrumentation  leads  through  magnetic  flux  lines  generated  by  the  DC 
motor,  (2)  motor  power  controller  noise  and  (3)  induced  alternating  currents 
through  the  lead  and  slipring  instrumentation  system  generating  fluctuating 
voltages.  The  voltage  data  used  in  the  data  reduction  program  was  obtained  by 
averaging  ten  successive  voltage  measurements  of  each  data  channel.  Repeat¬ 
ability  of  the  measurements  indicated  the  calculated  mean  temperature  was 
consistently  within  0.2®C  (0.36®F)  of  the  mean  temperature.  An  error  analysis 
of  the  data  reduction  equations  showed  that  approximately  3/4  of  the  estimated 
error  in  calculating  heat  transfer  coefficient  was  due  to  the  error  in  the 
temperature  measurement.  Estimates  in  the  error  in  calculating  heat  transfer 
coefficient  typically  varied  from  approximately  +2%  at  the  inlet  to  +10%  at 
the  exit  of  the  heat  transfer  model. 
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4.0  EXPERIMENTAL  PARAMETERS  AND  TEST  MATRIX 


The  present  study  of  heat  transfer  from  a  serpentine,  square-passage,  smooth 
wall  model  for  a  variety  of  stationary  and  rotating  flow  conditions  provided 
results  that  can  be  used  to  model  a  variety  of  geometric  locations  in  the 
internal  cooling  passages  of  a  gas  turbine  blade.  This  study  comprised  (i) 
experiments  originally  proposed  under  the  present  contract  (Tests  Nos.  1-13), 
(ii)  supplimentary  experiments  performed  under  United  Technologies  Corporation 
(P&W  and  UTRC)  independent  research  program  (Test  Nos.  101  -  117)  and  (iii) 
supplementary  experiments  performed  under  redirected  effort  on  the  present 
contract  (Test  Nos.  118-126).  The  test  conditions  for  these  experiments  are 
shown  on  Tables  I,  II,  and  III. 

A  dimensional  analysis  study  performed  at  UTRC  prior  to  the  onset  of  the 
present  study  (Suo,  1980),  similar  to  that  of  Guidez  (1988),  showed  that  the 
flow  patterns  and  hence  convective  heat  transfer  would  be  influenced  by  four 
nondimens ional  flow  parameters  and  several  geometric  parameters.  The 
nondimens ional  flow  parameters  are  as  follows; 

Reynold  number  -  pVd/fJ- 

Rotation  Number  -  Od/V 

Density  ratio  -  (P^,  -  Pw)/^b  ~  ('^w'^b^'^^w 

Buoyancy  Parameter  [  (Pb-Pv)/^bH  ^^R/V)  t  fl  d/V] 


For  flow  in  rotating  radial  coolant  passages,  Coriolis  forces,  represented  by 
the  nondimens ional  parameter, Hd/V,  and  the  nondimensional  streamwise  velocity 
gradients,  produce  secondary  flows  in  the  plane  perpendicular  to  the  radial 
direction.  These  secondary  flows  are  produced  by  the  viscous  force/Coriolis 
force  interaction.  Buoyancy  also  produces  secondary  flows  in  the  radial 
direction.  For  flow  in  rotating  radial  coolant  passages  with  walls  hotter  than 
the  bulk  fluid,  the  buoyancy  effects  always  tend  to  drive  the  heated  flow 
inward.  Thus  the  buoyancy  flow  direction  is  opposite  the  mean  velocity  direc¬ 
tion  for  flow  radially  outward  and  is  in  the  same  direction  for  flow  radially 
inward.  From  previous  studies,  both  the  Coriolis  and  buoyancy  forces  can  be 
expected  to  produce  significant  changes  to  the  coolant  passage  flow  field  and 
hence  heat  transfer.  Rotating  constant-temperature  flow  studies  by  Johnston  et 
al,  (1972)  have  shown  that  the  Coriolis  forces  can  deunpen  turbulent  fluctua¬ 
tions  and  lamlnarize  flow  in  portions  of  a  channel.  Combined  free  and  forced 
convection  studies  in  stationary  systems  have  shown  that  the  turbulent  shear 
structure  and  heat  transfer  is  significantly  altered  with  co-flowing  or 
counter-flowing  buoyancy  effects  (Eckert  et  al,  1953  and  Metals  and  Eckert, 
1964).  The  results  from  the  present  experimental  study  show  regions  where  the 
viscous,  Coriolis  or  buoyancy  forces  dominate  the  flow  field  and  regions  where 
the  interactions  between  the  forces  are  strong.  The  geometric  parameters  are 
as  follows: 


Axial  location 
Radial  location 
Flow  direction 
Passage  orientation 


-  X/d 

-  R/d 

-  inward,  outward 

-  a 
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TABLE  I 

TEST  CONDITIONS  FOR  ROTATING  HEAT  TRANSFER  EXPERIMENTS  WITH  SMOOTH  WALL  MODEL 

Contract  NAS3-23691 
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NOTES;  Re  =pVd/At  Gr/Re2  =  (AP>P  ) CfiR/V) (Qd/V) 

Ro  =nd/V  Gr  =  (AP/P)(ilR/V)(ncl/V)(pVd/jU)2 


Entrance  effects  (X/d)  are  expected  to  be  significant  in  each  of  the  three 
legs  as  the  flow  develops  downstream  of  the  inlet  and  downstream  of  the  turns. 
The  radial  location  (R/d)  of  the  model  was  varied  to  isolate  the  effects  of 
the  rotation  and  the  buoyancy  rotation  parameters.  The  passage  orientation 
was  varied  to  determine  effects  of  coolant  passage  orientation  at  the  leading 
and  trailing  regions  of  the  blade,  as  well  as  the  effects  of  the  secondary 
flow  interactions  caused  by  the  turns. 

The  variation  of  Reynolds  niunber,  rotation  number  and  inlet  wall  to  bulk 
temperature  difference  which  produced  the  parameteric  study  is  shown  in  Figure 
14.  The  standard  flow  condition  used  in  the  study  was  that  for  test  No.  4: 

Re  =  25,000,  Ro  =  fid/V  =  0.24,  ATin  =  44.4“C  (80“F),  R/d  =  49  and  O'  =  0.  The 
original  plan  was  to  vary  parameters  only  about  the  standard  flow  condition, 
varying  only  one  parameter  at  a  time.  However  at  the  completion  of  the  original 
test  plan,  it  became  apparent  that  the  heat  transfer  relationships  were  complex 
and  that  the  viscous,  Coriolis  and  buoyancy  forces  each  dominated  the  flow 
field  for  various  combinations  of  the  test  conditions. 

The  results  obtained  show  first  order  effects  for  the  following  parameters: 

1.  Streamwise  location  -  The  range  of  X/d  in  each  straight  passage  varies 
from  0  up  to  12.4  and  for  most  of  the  passage  length,  the  results  are 
in  the  developing  flow  region  for  constant  wall  temperature 
conditions.  However,  for  the  nonrotating  test  conditions  the  results 
approach  fully  developed  flow  levels  at  the  end  of  each  passage. 

2.  Reynolds  number  -  The  Reynolds  number  was  varied  from  12,500  to  50,000 
for  the  stationary  experiments  and  from  12,500  to  75,000  for  the 
rotating  experiments. 

3.  Rotation  Number  -  The  rotation  nioraber  Rd/V,  (the  inverse  of  the 
Rossby  number)  and  the  streamwise  velocity  gradients  are  the  primary 
nondimenslonal  factors  governing  secondary  flow  in  the  plane 
perpendicular  to  the  centerline  of  rotating  radial  ducts. 

4.  Density  ratio  -  The  density  ratio,  (  ”  ^w)'^  ^b* 

basic  nondlmehsional  parameters  obtained  from  several  previous  dimen¬ 
sional  analysis  of  flow  in  a  rotating  radial  duct.  The  product  of  the 
density  ratio  (^b'^w^/^b  ®  gravitational  pareuneter, 

( nd/V)2(R/d),  cause  secondary  flow  in  the  radial  direction.  For 
this  study  with  heated  walls  and  for  the  gas  turbine  blades,  the 
buoyancy  effect  is  always  radially  inward  whether  the  flow  direction 
is  radially  inward  or  outward.  Note:  (  ^b  “  ^w^^  ^b  =  ^^w  “  Tb)/T^, 
(AP//>)in  =  (AT/T)in. 

5.  Radius  ratio  -  The  ratio,  R/d,  was  also  obtained  from  several  dimen¬ 
sional  analysis  and  is  related  to  buoyancy  effects.  An  alternative 
nondimenslonal  parameter  could  be  used  as  the  fifth  basic  parameter, 
as  discussed  in  a  subsequent  paragraph.  However,  the  R/d  parameter  is 
readily  recognized  by  the  designer  and  researcher  alike. 
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ROTATION  NO.,  ROTATION  NO., 


a  =  0  (STD) 

R/d  =  49  (STD) 

NOTE:  STD.  BASELINE  FLOW  CONDITIONS 


REYNOLDS  NO.,  Re  REYNOLDS  NO.,  Re 


C)  AT|n  s:  66.700  (1200F)  D)  ATin*88.9®C  (160°F) 


REYNOLDS  NO.,  Re  REYNOLDS  NO.,  Re 


Figure  14.-  Test  Conditions  for  Parametric  Rotating  Heat  Transfer  Study 


6.  Passage  orientation  -  The  serpentine  model  was  constructed  such  that 
the  plane,  which  contains  the  centerlines  of  all  four  passages,  could 
be  rotated  about  a  radial  axis  through  the  geometric  centerline  of  all 
four  evenly-spaced  passages. For  Q!=  0“,  extensions  of  the  centerlines 
of  all  four  passages  would  pass  through  the  axis  of  rotation  as  shown 
in  Figure  9.  One  side  of  the  square  passage  becomes  the  leading  side, 
i.e.  nt.  For  a  =  45®,  each  test  section  passage  has  two  leading  and 
two  trailing  sides,  i.e.O’t* 

7.  Flow  Direction  -  The  direction  of  the  flow  causes  the  buoyancy, 
viscous  and  Coriolis  forces  to  interact  in  a  complex  manner.  Previous 
investigators  conducting  free  and  forced  convection  experiments  in 
stationary  tests  have  attributed  the  differences  in  heat  transfer 
between  flow  upward  and  downward  to  changes  in  the  turbulent  structure 
of  the  flow. 

A  total  of  thirty-nine  tests  were  conducted  with  the  smooth  wall  model  as 
shown  in  Tables  I,  II,  and  III  (Repeat  runs  were  also  obtained  under  the  same 
test  number).  In  order  to  make  the  presentation  of  the  principle  results  for 
this  program  tractable  and  discernible  to  the  reader,  the  heat  transfer  data 
is  presented  in  several  stages. 

1)  The  heat  transfer  results  are  presented  as  a  variation  of  the  basic 
parameters  about  a  standard  flow  condition  for  selected  streamwise 
locations. The  stationary  heat  transfer  results  is  presented  in  Section 
5.0  with  the  variation  of  Reynolds  number  and  density  ratio  about  the 
standard  stationary  flow  condition  (Test  1).  The  rotating  heat 
transfer  results  will  be  presented  in  Section  6  with  the  variation  of 
rotation  parameter,  density  ratio,  Reynolds  number,  model  radius  and 
model  orientation  angle  about  the  standard  rotating  flow  conditions 
(Test  4). 

2)  The  results  in  Section  7.0  is  presented  as  a  function  of  two  basic  or 
secondary  nondimensional  parameters  for  all  the  data  with  a  =  0®  . 

This  presentation  will  be  used  to  discern  the  complex  heat  transfer 
relationships  that  occurs  for  various  geometrical  locations  over  a 
wide  range  of  flow  conditions.  The  range  of  nondimensional  flow 
conditions  for  which  data  was  obtained  in  this  study  with  one  model 
geometry  encompasses  the  range  of  nondimensional  flow  conditions  in 
present  and  future  small  and  large  aircraft  gas  turbines. 

3)  Heat  transfer  correlations  based  on  the  results  of  this  experimental 
program  and  suitable  for  use  by  gas  turbine  designers  are  presented  in 
Section  8.0. 

4)  The  results  from  the  present  prograun  are  compared  in  Section  8.3  with 
results  from  previously  published  rotating  heat  transfer  experiments 
and  with  results  from  stationary  combined  free  and  forced  convection 
heat  transfer  experiments. 
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5.0  STATIONARY  HEAT  TRANSFER  RESULTS 


5.1  Baseline  Flow  Condition 

The  baseline  flow  condition  for  the  stationary  heat  transfer  experiments  is 
test  number  lA.  Test  lA  was  conducted  after  two-thirds  of  the  rotating  test 
program  was  conducted  and  after  the  heat  transfer  segments  had  "set"  due  to 
rotation.  Examination  of  the  smooth  wall  model  after  the  conclusion  of  all  the 
heat  transfer  experiments  showed  no  perceivable  changes  from  the  "before" 
condition.  Test  lA  was  used  as  the  basis  for  comparison  for  all  the  rotating 
tests  as  well  as  the  stationary  tests. 

The  variation  of  Nusselt  number,  Nu  =  hd/k,  with  streamwise  location  for  test 
lA  is  shown  in  Figure  15.  The  Nusselt  number  for  a  fully  developed  flow  in  a 
constant-wall-temperature,  square  duct  is  shown  for  comparison.  Note  that  the 
heat  transfer  rate  on  all  four  walls  decreases  by  a  factor  of  1.5  from  the 
first  to  the  third  heat  transfer  segment  in  each  straight  segment.  Nusselt 
numbers  for  most  of  the  third  segments  in  each  straight  passage  i.e..  Segments 
D,  I  &  N,  approach  the  Nusselt  number  for  fully  developed  flow  previously 
obtained  in  square  ducts. 

The  heat  transfer  increases  significantly  in  the  turn  sections,  as  expected. 
The  heat  transfer  on  the  leading  and  trailing  turn  surfaces  increases  a  factor 
of  approximately  two  greater  than  the  smooth  duct  value.  The  heat  transfer 
rate  is  approximately  10%  greater  in  the  second  half  of  the  turn  compared  to 
the  first  half.  The  heat  transfer  on  the  side  wall  surfaces  is  more  complex. 
The  first  half  of  side  walls  E,  J  and  P  are  continuations  of  the  straight 
section  side  walls.  The  second  half  of  side  walls  E,  J  and  P  and  the  first 
half  of  side  walls  F,  K  and  R  are  perpendicular  to  the  straight  segments.  The 
turbulence  level  of  the  flow  adjacent  to  these  surfaces  and  the  heat  transfer 
rate  are  both  expected  to  be  higher  than  the  fully  developed  flow.  The  second 
half  of  segments  F,  K  and  R  are  extensions  of  the  straight  passage  walls.  The 
heat  transfer  on  these  segments  is  high  but  somewhat  inconsistent  from  turn  to 
turn. 

The  conclusion  from  the  results  of  Test  lA  is  that  the  heat  transfer 
characteristic  for  this  smooth  wall  model  are  reasonably  well  behaved.  The 
heat  transfer  rates  on  each  of  the  four  walls  at  each  location  in  the  straight 
segments  are  generally  within  10  percent  of  the  average.  The  exception  is  for 
streamwise  location  D,  where  the  outside  sidewall  (test  surface  4)  increases, 
possibly  due  to  acceleration  effects  from  the  turn  and/or  separation  of  the 
flow  from  adjacent  surfaces. 

The  heat  transfer  results  for  the  first  two  straight  sections  are  represented 
in  Figure  16  as  the  ratio  of  local  Nusselt  number  to  the  Nusselt  number  for 
fully  developed  flow  in  a  square  duct  with  the  same  hydraulic  diameter  and  the 
same  level  Reynolds  Number.  This  heat  transfer  ratio,  Nu/Nu  will  be  used  in 
all  succeeding  comparisons.  Present  data  are  compared  with  previous  results 
from  entrance  region  heat  transfer  analysis  and  experiments  in  Bergles  and 
Webb  (1970).  The  results  from  the  present  experiment  lie  in  the  range  of 
previous  results.  Note  that  the  results  for  the  first  straight  passage  lie  in 
the  top  range  of  the  previous  results,  whereas,  those  from  the  second  straight 
passage  lie  in  the  center  to  lower  range  of  previous  results.  The  conclusion 
from  this  comparison  of  the  present  results  with  limited  results  from  previous 
experiments  is  that  the  present  data  generally  falls  in  the  expected  range. 
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Figure  15.-  Heat  Transfer  Results  for  Stationary  Baseline  Flow  Conditions 
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Figure  16.- 


Comoarlson  of  Heat  Transfer  in  Straight  Segments  for  Test 
Condition  Number  lA  With  Previous  Heat  Transfer  Results  for 
Entrance  Regions. 
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5.2  Effect  of  Density  Ratio 

The  effect  of  density  differences  between  the  wall  and  bulk  temperature  on  the 
heat  transfer  ratio  is  shown  in  Figure  17.  Results  are  presented  for  the 
baseline  flow  condition  (Test  lA)  with  a  wall  to  inlet  bulk  temperature 
difference  of  43.9°C  (79“F)  and  for  temperature  differences  of  22.2°C,  66.7°C, 
and  88.8®C  (40°F,  120‘’F  and  IbCF).  The  heat  transfer  ratios  are  essentially 
the  same  for  all  four  temperature  differences  at  each  heat  transfer  segment. 
The  fully  developed  Nusselt  number  used  for  comparison  was  calculated  at  each 
segment  using  the  local  film  temperature,  i.e.  the  average  of  the  bulk  and 
wall  temperature  for  the  determination  of  the  transport  properties  and  k) 
used  in  the  correlation.  The  major  conclusion  from  these  experiments  with 
various  wall  to  bulk  temperature  differences  is  that  the  use  of  the  film 
temperature  for  the  determination  of  transport  property  values  causes  the  heat 
transfer  ratio  results  to  be  independent  of  temperature  differences  for  this 
range.  This  conclusion  is  compatible  with  general  heat  transfer  practice  when 
the  heat  transfer  is  dominated  by  forced  convection. 

5.3  Effect  of  Reynolds  Number 

The  heat  transfer  ratio  distribution  for  three  Reynolds  numbers  is  shown  in 
Figure  18.  The  heat  transfer  ratio  at  each  streamwise  location  in  the  straight 
sections  is  generally  independent  of  Reynolds  number.  The  small  variations 
with  Reynolds  number  at  streamwise  locations  A  and  B  may  be  due  to  weak 
variations  in  the  entrance  region  length  with  Reynolds  number.  However,  the 
variations  are  of  the  same  order  as  the  experimental  uncertainty.  The  heat 
transfer  in  the  turns  shows  generally  higher  heat  transfer  ratios  at  the 
lowest  Reynolds  number  and  may  indicate  that  the  exponent  "N"  for  a  Nuoc  Re^ 
relationship  in  the  turns  is  less  than  the  value  0.8  used  for  the  Nu 
correlation.  The  conclusion  from  these  experiments  with  Reynolds  numbers 
equal,  above  and  below  the  stationary  baseline  flow  condition  is  that  the  heat 
transfer  ratio  for  stationary  conditions  is  generally  independent  of  Reynolds 
number  and  primarily  a  function  of  geometric  location. 

5.4  Concluding  Remarks 

The  heat  transfer  characteristics  of  the  smooth  wall,  serpentine  heat  transfer 
model  are  generally  well  behaved  and  in  agreement  with  previous  experiments 
and  heat  transfer  practice.  The  heat  transfer  results  in  each  straight  passage 
show  entrance  region  effects  and  are  approximately  equal  on  all  four  sides. 

The  heat  transfer  ratio  distribution  is  independent  of  bulk-to-wall 
temperature  differences  when  the  film  temperature  is  used  to  determine  the 
transport  properties.  The  heat  transfer  ratio  for  each  segment  is  a  weak 
function  of  Reynolds  number.  Most  of  the  differences  between  the  results  at 
various  Reynolds  numbers  was  within  the  expected  accuracy  for  these  heat 
transfer  experiments. 


29 
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Figure  17.-  Effect  of  Wall  Temperature  on  Heat  Transfer  Ratio  for  Nonrotating 
Flow  Conditions. 
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Figure  18.-  Effect  of  Reynolds  Number  on  Heat  Transfer  Ratio  for  Stationary 
Condition. 
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6.0  ROTATING  HEAT  TRANSFER  RESULTS 


6.1  Baseline  Flow  Condition 

The  baseline  flow  condition  for  the  rotating  experiments  is  test  number  4. 

Test  number  4  has  a  rotation  number,  fid/V,  equal  0.24,  a  temperature  ratio, 
(ziT/T)in  equal  to  0.13,  and  a  Reynolds  number  equal  to  25,000.  The  Reynolds 
number  and  the  density  ratio  are  the  same  as  those  for  the  stationary  base 
line  flow  condition.  The  rotation  number  and  the  ratio  of  buoyancy  forces  to 
viscous  forces  are  in  ranges  where  previous  constant-temperature,  rotating 
experiments  (Moon,  1964  and  Johnston,  et  al,  1972)  and  stationary  combined  free 
and  forced  convection  experiments  (Eckert,  et  al,  1953)  have  shown  strong  and 
moderate  effects,  respectively.  These  nondimens ional  test  flow  conditions  are 
also  in  the  central  region  of  the  nondimens ional  operating  flow  condition  range 
for  advanced  aircraft  gas  turbines.  In  addition,  the  rotating  baseline  flow 
condition  is  in  a  region  where  the  extrapolation  of  heat  transfer  correlations 
from  previous  rotating  heat  transfer  experiments  would  predict  large  decreases 
in  the  heat  transfer  coefficient  below  fully  developed  turbulent  rates  (Morris, 
1981).  Therefore,  at  least  modest  variations  in  the  heat  transfer  rate  are 
expected  to  occur  because  of  rotation. 

The  results  for  the  Rotating  Baseline  Flow  Condition  will  first  be  presented 
and  compared  with  the  results  for  the  Stationary  Baseline  Flow  Condition.  In 
the  following  paragraphs,  the  effects  of  rotation  and  buoyancy  will  be  related 
to  previous  rotating  flow  experiments  with  constant  density  fluids,  and  to 
stationary  heat  transfer  experiments  with  free  and  forced  convection  effects. 
The  author's  conjectures  will  be  identified  as  appropriate.  The  probable  cause 
and  effect  relationships  for  the  flow  and  heat  transfer  will  be  further 
discussed  as  the  experimental  results  for  variations  in  the  rotation  rate, 
density  ratio,  and  radius  ratio  are  presented. 

The  variations  of  heat  transfer  ratio  with  streeunwlse  location  are  shown  in 
Figure  19  for  the  rotating  and  stationary  base,  line  flow  conditions.  The  heat 
transfer  ratio  on  the  leading  and  trailing  segments  shows  the  largest 
variation  from  the  stationary  heat  transfer  values.  In  the  outward  straight 
passage  (streamwise  locations  A  to  D),  the  heat  transfer  ratio  decreases  to 
about  40%  of  the  stationary  value  on  the  leading  segment  at  streamwise 
location  C.  For  the  same  passage,  the  heat  transfer  ratio  increases  to  more 
than  2.3  times  the  stationary  value  on  the  trailing  segment  at  streamwise 
location  D.  On  both  sidewalls  in  the  first  outward  straight  passage,  the  heat 
transfer  ratio  increases  with  rotation,  compared  to  the  stationary  heat 
transfer  ratios. 

The  difference  in  heat  transfer  between  the  rotating  and  nonrotating  flow 
conditions  on  the  trailing  and  sidewall  surfaces  of  the  first  radially  outward 
flowing  passage  is  attributed  to  both  the  increasing  strength  of  the  secondary 
flow  cells  associated  with  the  Coriolis  force  and  the  buoyancy.  The  decrease 
in  heat  transfer  near  the  inlet  of  the  passage  on  the  leading  surfaces  is 
attributed  to  the  stabilizing  of  the  near-wall  flow,  as  observed  by  Johnston 
(1972).  The  subsequent  increase  in  heat  transfer  near  the  end  of  the  passage 
is  postulated  to  occur  when  the  secondary  flow  cells  become  more  developed  and 
interact  with  the  buoyant,  stabilized  near-wall  flow  on  the  leading  side  of 
the  passage.  Further  discussion  of  this  interaction  will  be  presented  in 
subsequent  sections. 
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Figure  19.-  Comparison  of  Heat  Transfer  Ratio  for  Rotating  and  Stationary 
Baseline  Flow  Conditions. 
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The  effect  of  rotation  on  heat  transfer  in  the  inward  straight  passage 
(streamwise  locations  G  to  I),  are  also  significant  but  less  dramatic.  The 
heat  transfer  ratio  increases  about  20%  on  the  leading  segments  compared  to 
the  stationary  results.  The  heat  transfer  on  the  trailing  segments  decreases 
from  5  to  30%  of  their  stationary  heat  transfer  rates.  The  heat  transfer  on 
the  side  walls  for  the  inward  straight  passage  is  generally  greater  with 
rotation  than  for  the  stationary  flow  condition.  The  exception  occurred  for 
segments  10  and  11  at  streamwise  locations  H  and  I  where  the  heat  transfer 
rates  are  5  percent  less  than  the  stationary  values. 

The  difference  in  heat  transfer  on  the  high  pressure  sides  of  the  coolant 
passage  (i.e.  trailing  surfaces  of  the  first,  outward  flowing  passage  and 
leading  surfaces  of  the  second,  inward  flowing  passage)  is  believed  to  be  a 
result  of  the  combined  effects  of  the  secondary  flow  induced  by  Coriolis 
forces  and  buoyancy  driven  secondary  flow.  When  a  counterflow  situation  exists 
(i.e.  buoyancy  driven  flow  in  a  direction  opposite  the  mean  flow  direction  as 
occurs  for  outward  flow  and  heated  walls),  the  combined  effects  of  buoyant  and 
Coriolis-driven  secondary  flows  cause  significant  increases  in  the  heat 
transfer  coefficients.  When  a  parallel  flow  situation  exists  (i.e.  buoyancy 
driven  flow  in  the  same  direction  as  the  mean  flow  direction),  the  combined 
effects  of  rotation  and  buoyancy  appear  to  counteract  each  other  and  only  a 
modest  (10  percent)  increase  in  heat  transfer  occurs,  compared  to  the 
stationary  results.  Futher  discussion  of  these  effects  will  be  presented  in 
subsequent  sections. 

The  heat  transfer  in  the  turn  regions  is  also  significantly  increased  and 
decreased  by  rotation.  Note  that  the  heat  transfer  ratios  for  the  first 
outside  turn  (streamwise  locations  E  and  F)  Increases  by  50%  on  all  segments 
compared  to  the  nonrotating  values.  The  heat  transfer  rates  for  the  second 
inside  turn  (streamwise  locations  J  and  K)  decreases  by  an  amount  equal  to  15 
to  50%  of  the  stationary  value.  Because  the  Increased  heat  transfer  on  the 
outside  turn  and  decreased  heat  transfer  on  the  Inside  turn  occurs  on  both  the 
leading  and  trailing  surfaces,  these  changes  aije  more  likely  to  be  caused  by 
buoyancy  effects  rather  than  by  the  conservation  of  vorticity  through  the 
turn.  Note,  however,  that  the  heat  transfer  for  the  outside  turn  at  streamwise 
location  F  is  considerably  greater  for  the  leading  segment  than  the  trailing 
segment.  These  asymmetric  differences  in  the  heat  transfer  ratio  are  likely  to 
be  caused  by  the  convection  of  the  secondary  flow  (produced  in  the  straight 
passages  by  Coriolis  forces)  around  the  180*  turn. 

The  conclusion  from  this  comparison  of  the  heat  transfer  ratios  for  the 
rotating  and  stationary  flow  conditions  is  that  the  effects  of  rotation  on 
heat  transfer  in  a  radial  passage  are  significant.  The  heat  transfer  rate 
decreased  to  less  than  half  the  stationary  value  on  a  leading  wall  segment, 
increased  by  a  factor  of  two  on  a  trailing  wall  segment.  Increased  by  as  much 
as  a  factor  of  two  on  the  outside  turn  and  decreased  by  as  much  as  50%  in  the 
inside  turn.  The  presentation  of  results  at  other  rotating  flow  conditions, 
obtained  by  varying  one  of  the  four  basic  dimensionless  parameters,  will 
isolate  the  effects  of  each  parameter.  The  effects  of  the  rotation  numbers, 
density  ratio,  Reynolds  niimbers  and  radius  ratio  will  be  discussed  in  the 
aforementioned  order  which  is  also  the  order  of  decreasing  importance  to  heat 
transfer. 
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6.2  Effect  of  Rotation  Number 

The  effect  of  rotation  number,  Od/V,  on  heat  transfer  ratio  is  shown  in  Figure 
20.  All  other  flow  conditions  i.e.,/^T,  R/d,  Re  and  a,  are  held  constant  at 
their  baseline  values,  illustrative  results  are  presented  for  5  rotation  rates 
including  the  stationary  and  rotating  baseline  flow  conditions.  Additional 
data  are  available  at  15,  137  and  380  rpm.  The  effects  of  rotation  on  the  heat 
transfer  ratio  are  complex  for  some  geometric  regions  and  are  monotonic  in 
others . 


I 

6.2.1  High  Pressure  Surfaces 

The  most  straightforward  effect  of  rotation  is  observed  on  the  trailing 
segments  of  the  first  outward  straight  passage  (streamwise  locations  A  to  D). 

On  these  segments,  the  heat  transfer  ratio  increases  monotonically  with 
rotation  at  all  streamwise  locations.  Note  that  the  largest  relative  increases 
occur  at  streamwise  locations  C  and  D  for  a  rotation  number  change  from  0  to 
0.12.  The  largest  absolute  increase  occurs  at  streamwise  location  B  and  C  for 
a  rotation  number  change  from  0.24  to  0.35.  Note  that  the  heat  transfer  ratio 
increases  by  more  than  a  factor  of  3.5  at  streamwise  location  C  for  rotation 
number  of  0.48  compared  to  the  zero  rotation  value.  Smaller  increases  in  the 
heat  transfer  ratio  on  the  leading  segments  occurred  for  the  inward  flow 
straight  passage  (streamwise  locations  G,  H  and  I). 

The  effects  on  heat  transfer  due  to  Coriolis  generated  secondary  flows  might 
be  expected  to  be  approximately  the  same  for  the  trailing  segments  of  the 
first  passage  and  the  leading  segments  of  the  second  passage.  The  differences 
in  heat  transfer  between  the  outward  and  inward  flowing  passages  are  therefore 
attributed  to  the  different  effects  of  buoyancy  in  the  counter-flowing  first 
passage  (radially  outward  flow)  and  the  co-flowing  second  passage  (radially 
inward  flow) . 

The  lesser  increases  in  the  heat  transfer  ratio  on  the  high  pressure  side  of 
the  second  passage  are  attributed  to  the  counteracting  effects  of  rotation  and 
parallel-flow,  combined-free-and-forced  convection  in  the  passage.  Previous 
authors  have  attributed  the  decreased  heat  transfer  in  stationary  parallel- 
flow  free-and-forced-convection  experiments  to  a  reduction  in  the  generation 
of  near-wall  turbulence.  In  the  first  passage,  the  near-wall  buoyancy  driven 
flow  was  inward  toward  the  axis  of  rotation  and  the  coolant  flow  was  outward. 
This  counter  flow  situation  generated  additional  near-wail  turbulence  due  to 
the  strong  shear  gradient.  This  increase  in  the  shear  forces  combined  with  the 
cross  stream  secondary  flows  generated  by  Coriolis  forces  cause  large  increases 
in  heat  transfer  in  the  first  passage.  However,  when  the  flow  and  the  buoyancy 
driven  near-wall  flows  are  coincident,  as  in  the  second  passage,  the  generation 
of  near-wall  turbulence  may  be  decreased  because  of  the  changes  in  the  velocity 
profile.  Therefore,  the  combined  effects  of  the  buoyant  and  the  cross  stream 
secondary  flows  in  the  second  passage  on  the  heat  transfer  are  less.  The 
magnitude  of  the  buoyancy  effect  on  the  heat  transfer  is  unclear  in  that  the 
buoyancy  effect  in  the  second  passage  may  be  zero  (which  implies  a  modest 
Coriolis  dominated  heat  transfer  increase)  or  negative  (which  implies  a  larger 
Coriolis  dominated  heat  transfer  increase  which  is  offset  by  a  reduction  due 
to  buoyancy).  The  cause  and  effect  relationships  presented  in  this  paragraph 
are  the  author's  speculations,  and  will  require  further  experimentation  or 
calculations  by  direct  numerical  simulation  in  order  to  be  substantiated  or 
rejected. 
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Figure  20.-  Effects  of  Rotation  Rate  on  Heat  Transfer  Ratio. 
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6.2.2  Low  Pressure  Surfaces 


In  contrast  to  the  continual  increase  in  heat  transfer  with  increasing 
rotation  number  on  the  trailing  side  of  the  first  passage,  the  heat  transfer 
ratio  decreases  with  increasing  rotation  number  on  the  leading  side  of  the 
passage  near  the  inlet.  For  all  of  the  remaining  locations  on  the  leading  side 
of  the  passage,  the  heat  transfer  ratio  decreases  and  then  Increases  again 
with  increasing  rotation  number.  Heat  transfer  from  the  trailing,  low  pressure 
surfaces  of  the  second  passage  also  had  large  decreases  in  heat  transfer.  Heat 
transfer  in  the  second  passage  decreased  to  almost  60%  of  the  stationary  heat 
transfer  levels  compared  to  40%  in  the  first  passage.  In  the  second  passage, 
the  heat  transfer  decreased  and  then  subsequently  increased  again  as  the 
rotation  rate  was  increased. 

The  decreases  in  the  heat  transfer  ratio  are  attributed,  for  the  most  part,  to 
the  Coriolis  generated  cross  stream  flow  patterns  as  well  as  the  stabilization 
of  the  near-wall  flow  on  the  leading  side  of  the  first  passage  (Johnston  et 
al.,  1972).  The  cross  stream  flows  cause  heated,  near-wall  fluid  from  the 
trailing  and  sidewall  surfaces  to  accumulate  near  the  leading  side  of  the 
coolant  passage  resulting  in  reduced  heat  transfer.  In  addition,  the  rotation 
stabilizes  the  shear  layers  along  this  wall  and  further  reduces  the  turbulent 
transport  of  heat.  The  Increase  in  the  heat  transfer  ratio  in  the  latter  half 
of  the  coolant  passage  for  the  larger  rotation  numbers  is  attributed  to 
increases  in  the  buoyancy  effects.  These  could  include  the  formation  of  radial 
recirculation  cells.  Similar  effects  of  rotation  are  noted  for  the  low 
pressure  surfaces  in  both  the  first  and  second  passages.  Irrespective  of  flow 
direction.  These  results  suggest  that  the  heat  transfer  on  low  pressure 
surfaces  is  dominated  by  Coriolis  generated  cross  stream  flows  which  cause  a 
stabilization  of  the  near-wall  flows  and  that  the  heat  transfer  on  the  high 
pressure  surfaces  is  affected  by  a  combination  of  Coriolis  and  buoyant 
effects.  Therefore,  it  can  be  expected  that  the  correlations  of  local  heat 
transfer  data  may  be  substantially  different  depending  on  local  flow 
conditions  (i.e.  due  to  differing  near-wall  shear  gradients). 

6.2.3  Side  Walls 


The  effects  of  rotation  on  the  heat  transfer  from  the  side  walls  in  the 
straight  passage  are  less  than  the  leading  and  trailing  segments.  The  general 
effect  is  that  the  side  wall  heat  transfer  increases  with  rotation.  The 
increases  in  heat  transfer  ratio  on  the  side  wall  segments  (1-4  and  19-22)  in 
the  first  passage  with  outward  flow  are  approximately  one-third  the  Increases 
on  the  trailing  segments  (49-52).  The  Increases  on  the  side  wall  segments 
(9-11  and  23-25)  in  the  second  straight  section  with  inward  flow  are  equal  or 
less  than  the  Increases  on  the  leading  segments  (39-41). 

6.2.4  Turns 


The  effect  of  rotation  number  on  heat  transfer  from  the  outside  turn 
(streamwise  locations  E  and  F)  is  monotonic  for  almost  all  points;  the  heat 
transfer  ratio  Increases  with  rotation  rate.  The  effect  of  rotation  number  on 
heat  transfer  ratio  from  the  inside  turn  (streamwise  locations  J  and  K)  is 
more  complex.  For  the  inside  turn,  the  heat  transfer  decreases  on  the  trailing 
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segments  (58,  59)  and  side  walls  (26-29)  as  the  rotation  number  is  increased 
from  0  to  0.24  and  then  increases  slightly  as  Ro  is  increased  to  0.48.  The 
heat  transfer  on  the  leading  segments  (42,  43)  decreases  from  0  to  20%  as  Ro 
is  increased  to  0.24  and  then  increases  3,0%  as  Ro  is  increased  from  0.24  to 
0.48. 


Although  the  rotation  number  is  a  basic  dimensionless  parameter  governing  the 
flow  in  rotating  passages,  the  variations  of  the  rotation  number  in  these  heat 
transfer  experiments  also  causes  buoyancy  effects.  Therefore,  no  conclusions 
regarding  the  effects  of  rotation  number  alone  can  be  deduced  without  the 
analysis  of  further  experiments  where  the  buoyancy  effects  are  also  fixed.  The 
effects  of  rotation  number  are  isolated  in  the  analysis  of  results  from  a  more 
complete  set  of  flow  conditions  in  Section  7. 

6.3  Effect  of  Density  Ratio 

6.3.1  Straight  Sections 

The  effects  of  varying  the  density  ratio  (wall  temperature)  on  the  heat 
transfer  ratio  are  shown  in  Figure  21.  All  flow  conditions  were  held  at  the 
rotating  baseline  flow  condition  except  the  inlet  density  ratio  (AP/P  )ij^ 
which  was  varied  from  0.07  to  0.22.  These  sets  of  data  were  obtained  at  a 
constant  rotation  number  and  therefore  conclusions  can  be  obtained  regarding 
the  effects  of  buoyancy  for  flow  conditions  near  the  rotating  baseline  flow 
conditions . 

For  the  trailing  and  side  wall  segments  in  both  the  outward  and  inward  flow 
straight  passages,  increasing  the  density  ratio  (and  hence  buoyancy)  increases 
the  heat  transfer  ratio.  This  increase  ranges  from  10  to  50%,  depending  upon 
location. 

Increasing  the  inlet  density  ratio  (i.e.  the  wall-to-coolant  temperature 
difference)  from  0.07  to  0.22  causes  the  heat  transfer  ratio  in  the  first 
passage  to  increase  on  the  trailing  surfaces  by  as  much  as  50%  and  on  the 
leading  surfaces  by  as  much  as  100%.  The  exception  to  the  general  increase  in 
heat  transfer  with  increasing  density  ratio  occurred  near  the  inlet  of  the 
first  passage  on  the  leading  side,  where  the  heat  transfer  ratio  is  observed 
to  decrease  slightly. 

Heat  transfer  in  the  second,  inward  flowing  passage  also  increases  with 
increasing  density  ratio.  In  general,  the  increases  in  heat  transfer  in  the 
second  passage  were  approximately  half  of  those  in  the  first  passage  (on  the 
order  10  to  50%  compared  to  maximum  relative  increase  of  100%  in  the  first 
passage) . 

The  differences  in  heat  transfer  behavior  due  to  changes  in  the  density  ratio 
between  the  first  and  second  passages  are  attributed  to  the  differing 
mechanisms  of  Coriolis  and  buoyancy  interaction.  If  the  effect  of  Coriolis 
generated  secondary  flow  on  heat  transfer  is  similar  (regardless  of  flow 
direction)  and  the  effect  of  increasing  density  ratio  for  fixed  rotation 
number  generally  causes  heat  transfer  to  increase,  then  the  interaction  of  the 
two  effects  is  significant  and  also  counteracting.  The  counteraction  of  the 
two  effects  was  evident  in  the  relatively  small  increases  in  heat  transfer  on 
the  high  pressure  side  of  the  second  passage.  The  reason  for  this  behavior  is 
not  known  at  this  time. 
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Figure  21.-  Effects  of  Density  (Temperature)  Ratio  on  Heat  Transfer  Ratio. 
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6.3.2  Turns 


The  variation  of  the  heat  transfer  ratio  with  increasing  density  ratio  was 
montonic  in  both  the  outside  and  inside  turns.  Increasing  the  density  ratio 
(wall-to-bulk  temperature  difference)  causes  the  heat  transfer  from  all 
surfaces  on  the  outside  turn  (locations  E  and  F)  walls  to  Increase.  However, 
increasing  the  density  ratio  causes  the  heat  transfer  from  all  surfaces  on  the 
inside  turn  (locations  J  and  K)  to  decrease.  Because  the  heat  transfer  on  the 
leading  and  trailing  sides  of  the  turns  (for  both  the  first  and  second  turns) 
is  not  equal,  and  because  Increasing  the  density  ratio  does  not  significantly 
reduce  the  differences,  it  is  plausible  that  the  transport  of  the  secondary 
flow  vorticity  from  the  upstream  straight  sections  has  an  important  effect  on 
heat  transfer  ratio  near  the  rotating  baseline  flow  condition. 


The  conclusion  from  this  discussion  is  that  increasing  the  density  ratio  (wall 
to  bulk  temperature  difference)  generally  causes  the  heat  transfer  in  both  the 
inward  and  outward  flow  straight  passages  to  increase.  Heat  transfer  in  the 
turns  increases  on  the  outside  turn  and  decreases  on  the  inside  turn  as  the 
density  ratio  increases. 

6.4  Effect  of  Reynolds  Number 

The  effect  of  varying  the  Reynolds  number  above  and  below  the  baseline  value 
while  maintaining  all  other  basic  dimensionless  parameters  at  the  rotating 
baseline  flow  conditions  value  is  shown  in  Figure  22.  In  the  first  passage 
(location  A  to  D),  the  heat  transfer  ratio  increases  modestly  (10%)  on  most 
segments  for  both  an  increase  and  decrease  in  Reynolds  number.  In  the  second 
passage  (location  G,  H  and  I),  the  increase  was  somewhat  greater  (up  to  25%). 

The  variation  of  heat  transfer  ratio  with  Reynolds  number  at  the  rotating  flow 
conditions  is  considerably  greater  than  observed  for  the  same  variation  at 
their  stationary  flow  conditions.  Although  the  basic  dimensionless  parameters 
are  constant  and  thus  the  ratio  of  Grashoff/Reynolds^,  a  parameter  judged  to 
influence  the  relative  effects  of  free  and  forced  convention,  is  constant  for 
these  experiments,  the  absolute  value  of  the  Grashoff  number  varies  by  a 
factor  of  16  for  a  variation  of  the  Reynolds  number  by  a  factor  of  4  (from 
12,500  to  50,000).  The  conclusion  from  these  tests  are  that  the  Reynolds 
number  effects  are  generally  contained  in  the  Nu^^jOc  Re®*®  parameter  used 
in  the  denominator  of  the  heat  transfer  ratio  and  that  the  use  of  the  heat 
transfer  ratio  is  an  acceptable  method  of  normalizing  the  test  results  for 
these  rotating  heat  transfer  experiments.  The  effects  of  Reynolds  number 
variation  are  also  discussed  in  Section  8.0  -  Correlation  of  Heat  Transfer 
Results. 

6.5  Effect  of  Model  Radius 

The  effects  of  buoyancy  are  coupled  as  the  product  of  three  of  the  four  basic 
flow  parameters  into  a  combined  buoyancy  parameter:  (S2d/V)2(R/d) (/AP/P ) .  In 
order  to  isolate  the  effects  of  R/d  and  obtaiij  heat  transfer  data  at  several 
conditions  with  the  same  product  of  (R/d)(AP/P),  the  model  radius  was  decreased 
to  about  two-thirds  of  its  baseline  value.  The  effects  of  this  decrease  in 
model  radius  on  the  heat  transfer  ratio  are  shown  in  Figure  23.  The  overall 
impression  is  that  decreasing  the  radius  did  not  significantly  change  the  heat 
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Figure  22.-  Effects  of  Reynolds  Number  on  Heat  Transfer  Ratio. 
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Figure  23.-  Effects  of  Model  Radius  on  Heat  Transfer  Ratio. 
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transfer  distribution.  However,  there  are  several  effects  which  are  consistent 
with  previously  described  results.  The  heat  transfer  ratios,  Nu/Nu^^,  on  both 
the  leading  and  trailing  segments  at  location  C  and  D  are  slightly  less  (0.1 
to  0.2  Nu/NUqq  )  than  the  rotating  baseline  values.  These  results  are  similar 
to  those  discussed  in  section  6.3  where  density  ratio  was  decreased  from  the 
baseline  value  of  0.13  to  0.07.  Because  the  effects  of  buoyancy  are  coupled  as 
the  product  of  two  flow  parameters  and  one  geometric  parameter  into  a  combined 
buoyancy  parameter  (AP/P ) (R/d) (nd/V)2 ,  varying  the  density  ratio  or  the 
radius  ratio  by  similar  amounts  should  cause  similar  variations  in  the  heat 
transfer  distributions. 

6.6  Effect  of  Model  Angular  Orientation 

The  baseline  rotating  flow  condition  has  the  test  model  oriented  such  that  the 
centerlines  of  all  four  straight  passages  pass  through  the  axis  of  rotation. 
However,  coolant  passages  are  typically  oriented  at  various  angles  to  the  axis 
of  rotation,  and  therefore,  it  is  desirable  to  determine  these  effects.  The 
model  was  rotated  45“  on  the  support  arm  as  shown  in  Figure  9.  The  result  is 
that  segments  33  to  36  and  19  to  22  are  co-leading  segments  in  the  first 
passage  and  segments  1  to  4  and  49  to  52  are  the  co-trailing  segments.  The 
effect  of  model  orientation  on  the  heat  transfer  ratio  is  shown  in  Figure  24. 
All  the  flow  conditions  for  test  no.  12  are  nearly  identical  to  those  for  test 
no.  4  (the  rotating  baseline  flow  conditions),  except  for  the  model 
orientation  angle  «  .  Note  that  for  the  first  outward  passage  of  test  no.  12, 
the  heat  transfer  from  the  co-trailing  segments  49-52  and  segments  1-4  for 
streamwise  locations  A,  B,  C  and  D  are  approximately  equal.  Likewise,  the  heat 
transfer  from  the  co-leading  segments  33-36  and  segments  19-22  are 
approximately  equal.  A  similar  trend  is  shown  in  the  first  inward  passage. 
Thus,  the  heat  transfer  is  symmetric  about  a  diagonal  in  the  direction  of 
rotation  across  the  passage.  The  heat  transfer  on  both  the  co-leading  and 
co-trailing  segments  (  a  =  45“ )  is  greater  than  or  equal  the  heat  transfer 
from  the  leading  and  trailing  segments  (a  =  0“)  for  the  baseline  flow 
condition,  with  the  exception  of  Segments  19-29  which  show  lower  heat  transfer 
for  a  =  45“  than  fora  =  0“ . 

The  effects  of  model  orientation  on  the  heat  transfer  ratio  shown  in  Figure  24 
are  also  shown  in  Figures  25  and  26  as  the  variation  of  heat  transfer  around 
the  test  section  at  selected  streamwise  locations.  The  results  from  all  four 
streamwise  locations  in  the  first  outward  sections  are  presented  in  Figure  25; 
the  results  for  two  streamwise  locations  in  all  three  legs  are  compared  in 
Figure  26.  The  symmetry  of  the  heat  transfer  ratio  for  fld/V  =  0.24  noted  for 
the  first  leg  is  apparent  in  these  presentations.  Small  asymmetries  occur  in 
the  second  and  third  leg  with  Hd/V  =  0.24. 

The  combined  effects  of  model  orientation  and  rotation  number  on  the  heat 
transfer  ratio  in  the  first  leg  are  also  shown  in  Figure  25  (ForHd/V  =  0, 
heat  transfer  data  obtained  for  a=  0“  is  also  replotted  at  the  same  segment 
location  for  a=  45“ ).  Note  the  change  of  shape  of  the  heat  transfer  ratio 
distributions  as  the  flow  progresses  from  the  guard  (streamwise  location  A)  to 
the  third  section  (streamwise  location  D)  in  Figure  25.  A  key  result  is  that 
the  decrease  in  heat  transfer  ratio  on  the  leading  surfaces  due  to  rotation  is 
less  for  a=  45“  than  for  Q!=  0“.  The  minimum  heat  transfer  ratio  in  the  first 
leg  for  a  =  45“  was  approximately  0.9,  whereas,  the  minimum  value  of  the  heat 
transfer  ratio  for  0=0“  was  0.45. 
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Figure  24.-  Effects  of  Model  Orientation  on  Heat  Transfer  Ratio. 
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The  combined  effects  of  model  orientation  and  rotation  number  on  the  heat 
transfer  ratio  in  all  three  legs  are  shown  in  Figure  26.  The  heat  transfer 
ratio  distributions  for  the  second  and  third  legs  are  somewhat  less  symmetrical 
than  the  distributions  for  the  first  leg.  The  heat  transfer  ratio  distribution 
for  streamwise  location  I  appears  anomalous.  Note  that  at  this  location,  the 
heat  transfer  ratios  on  the  co-leading  surfaces  for  a  =  45“  are  greater  than 
the  ratios  for  on  the  leading  surface  for  a;  =  0“ .  This  phenomena  occurred  only 
at  this  streamwise  location  where  the  flow  was  inward. 

The  results  from  the  experiments  with  a  =  45“  were  well  behaved  in  the  first 
leg  of  the  model  and  more  complex  in  the  second  and  third  leg.  The  heat 
transfer  ratio  distributions  in  the  second  and  third  leg  are  less  well 
behaved,  probably  due  to  the  more  complex  interactions  as  the  flow  passes 
through  the  turns  and  becomes  reestablished  in  the  straight  sections.  Although 
the  amount  of  data  available  from  the  a  =  45“  was  limited,  the  results  do  show 
the  general  effects  of  model  orientation. 

6.7  Conclusions 

The  heat  transfer  characteristics  of  the  smooth  wall  serpentine  heat  transfer 
model  with  rotation  are  considerably  more  complex  than  without  rotation.  The 
results  showed  large  effects  of  rotation  number  on  the  heat  transfer  ratio  for 
the  leading  and  trailing  surfaces.  In  the  first  leg  of  the  model,  the  heat 
transfer  ratio  decreased  to  as  much  as  40%  of  the  stationary  value  on  the 
leading  surfaces  and  increased  to  as  much  as  300%  of  the  stationary  value  on 
the  trailing  surfaces.  This  resulted  in  a  factor  of  7.5  difference  in  the  heat 
transfer  coefficients  on  opposite  sides  of  the  coolant  passage.  The  effects  of 
density  ratio  were  also  significant  although  less  than  the  effects  of  rotation 
number.  In  all  of  the  straight  sections,  increasing  the  density  ratio  caused 
the  heat  transfer  ratio  to  Increase.  (The  sole  exception  was  for  the  guard 
heater  on  the  leading  segment  in  the  first  leg.)  The  increase  in  the  heat 
transfer  ratio  was  as  much  as  100%  for  an  increase  in  the  inlet  wall-to-bulk 
temperature  difference  from  22.2“C  to  88.9“C  (40“F  to  160“F).  This  corresponds 
to  an  increase  in  the  inlet  density  ratio  from  0.07  to  0.22.  The  effects  of 
Reynolds  number  were  less  than  those  of  rotation  number  and  density  ratio.  The 
effects  of  model  radius  location  were  also  moderate  and  will  be  discussed 
further  in  the  next  section.  The  effect  of  model  orientation  was  also  complex. 
One  beneficial  result  of  model  orientation  (0;=  45“)  was  that  the  very  low 
heat  transfer  ratios  measured  on  the  leading  surface  in  the  first  leg  were 
ameliorated. 


47 


7.0  DISCUSSION  OF  COMBINED  EFFECTS  ON  HEAT  TRANSFER 


A  dimensional  analysis  of  flow  in  radial  rotating  ducts,  conducted  prior  to 
the  onset  of  this  study  (Appendix  10.1),  showed  that  two  combinations  of 
dimensional  parameters  cause  changes  in  a  flow  field  within  a  rotating  heated 
duct,  compared  to  the  stationary  duct.  There  are: 

(O  d/V)(l  +  CAP/P) )/(2  T  x'v)  (7.1) 

Coriolis  forces  causing  secondary  flow  in  the  plane  perpendicular  to  the  flow 
direction,  and 

(H  d/V)(AP/P)(nR/V)r  (7.2) 

buoyancy  forces  causing  secondary  flow  in  the  radial  direction. 

In  the  limit  as  A  P/P-0,  only  the  Coriolis  forces  are  present.  For  this  case, 
the  Reynolds  number  (PVd/P)  and  the  rotation  number  (Hd/V)  should  correlate 
results  for  a  given  geometry.  As  the  parameter  (AP/P)(nR/V)  approaches  the 
value  of  1,  Coriolis  forces  and  buoyancy  forces  are  of  equal  magnitude  and 
will  begin  to  interact  in  a  complex  manner.  When  the  parameter  (AP/P ) (flR/V) 
becomes  order  of  10,  buoyancy  forces  will  tend  to  dominate  the  flow  field. 

Early  rotating  heat  transfer  studies,  e.g.  Morris  et  al,  (1979)  used  relation¬ 
ships  of  the  type;  Nu  =  A  (Gr/Re2)Bpr%o*  to  correlate  their  results. 

(The  parameter  Gr/Re^  is  equivalent  to  (nd/V)(fi  R /V) (AP/P)  or  (nd/V)2 
(R/d)(AP/P)i  the  aforementioned  buoyancy  parameter.  From  the  previous 
paragraph,  it  can  be  discerned  that  Coriolis  and  buoyancy  forces  are  closely 
coupled  through  the  presence  of  the  rotation  number.  In  order  to  understand 
the  complex  cause/effect  relationships  between  forced  convection,  Coriolis  and 
buoyancy  forces,  the  results  from  the  experiments  will  first  be  presented  as 
functions  of  the  basic  parameters,  fid/V  and  (AP/P)in»  then  be  presented 
as  functions  of  the  buoyancy  parameter  with  the  rotation  pareuneter  noted.  In 
this  manner,  the  flow  regimes  where  viscous,  Coriolis  and  buoyancy  forces 
dominate  and  interact  will  be  identified. 

7.1  Effects  of  Density  Ratio 

The  variations  of  heat  transfer  ratio  with  inlet  density  ratio  (AP/P)in 
the  test  surfaces  on  the  leading  and  trailing  sides  of  the  first  leg  are  shown 
in  Figures  27  and  28  for  Re  =  25,000  and  selected  rotation  numbers,  fld/V. 
Curves  have  been  drawn  through  sets  of  data  at  the  same  rotation  number  and 
extrapolated  to  (AP/P)in  =  0  to  determine  the  effects  of  Coriolis  forces 
without  buoyancy  effects.  The  extrapolations  for  the  trailing  side  of  the 
model,  (Figure  28  -  Heaters  50,  51  and  52)  are  fairly  linear.  Note  that,  for 
fld/V  =  0,  the  heat  transfer  ratio  is  independent  of  (AP/P)in  expected 
when  Nu/NUj^j  is  determined  using  film  properties  (see  Section  5.2). 
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Determination  of  the  variation  of  Nu/NUo^with  (^P/P)in  near  0  is  more 
difficult  for  the  data  from  the  leading  test  surface  (Figure  27  -  Heaters  34, 
35  and  36)in  the  first  leg.  Direct  extrapolation  of  the  data  forOd/V  >  0.18 
(SymbolsD  and  -  Heater  36)  and  (AP/P)in  <  0.1  with  large  negative 
curvatures  would  lead  to  estimated  values  of  Nu/NUo^at  (AP/P)in  =  0  which 
are  less  than  0.35,  the  minimum  value  measured  thus  far.  Therefore,  the 
"dashed"  curve  leading  for  the  extrapolated  values  of  Nu/NUoohave  been 
"forced"  to  values  greater  than,  or  equal  to  0.35  at  (  AP/P)in  =  0.  It  is 
believed  that  the  secondary  flow  produced  by  the  Coriolis  forces  will  cause 
the  heat  transfer  rates  to  reach  a  minimum  value,  i.e.  Nu/NUg^=  0.35,  and 
then  remain  constant  or  possibly  increase. 

7.2  Effect  of  Rotation  Number 

Additional  data  from  parametric  variations  of  density  ratio  and  rotation 
parameter  were  necessary  to  isolate  the  effects  of  rotation  and  buoyancy.  The 
inlet  density  ratio  was  varied  from  0.07  to  0.22  for  selected  rotation 
numbers.  Heat  tansfer  results  from  these  experiments  were  plotted  vs.  inlet 
density  ratio  with  rotation  number  as  a  secondary  variable  (section  7.1).  The 
distributions  of  heat  transfer  ratio  with  density  ratio  were  extrapolated  for 
each  value  of  the  rotation  number  to  obtain  a  value  of  the  heat  transfer  ratio 
for  a  density  ratio  of  0.0  (i.e.  limit  as  AT  approaches  0.0).  The  heat 
transfer  results  obtained  from  the  experiments  plus  the  extrapolated  values 
for  a  density  ratio  of  0.0  (dashed  lines)  are  presented  in  Figures  29  and  30 
as  the  variation  of  the  heat  transfer  ratio  with  the  rotation  niimber  with  the 
inlet  density  ratio  as  the  secondary  variable  for  three  streamwise  locations 
for  the  first  and  second  passage. 


7.2.1  High  Pressure  Surfaces 

Heat  transfer  results  from  the  high  pressure  side  of  the  first  and  second 
passages  are  shown  in  Figures  29  and  30  for  ranges  of  rotation  number  and 
inlet  density  ratio.  Note  that  there  is  no  effect  of  density  ratio  on  the  heat 
transfer  ratio  for  a  rotation  number  of  0  when  film  properties  are  used  for 
the  dimensionaless  heat  transfer  and  flow  parameters.  Increasing  the  rotation 
number  causes  local  increases  in  the  heat  transfer  in  the  first  passages  by  as 
much  as  3.5  compared  to  the  heat  transfer  for  a  rotation  niimber  of  0.  Whereas 
the  heat  transfer  ratios  for  the  high  pressure  surfaces  increase  sharply  with 
increases  in  either  the  density  ratio  or  the  rotation  niamber,  with  one 
exception,  heat  transfer  in  the  second  passage  is  relatively  unaffected  by 
variations  of  either  parameter.  The  exception  being  near  the  inlet  of  the 
second  passage,  just  downstream  of  the  first  turn.  At  this  location,  the  heat 
transfer  increases  slightly  with  increases  in  the  rotation  parameter  and  the 
density  ratio.  However,  for  larger  X/d  in  the  second  passage,  the  effect  on 
the  heat  transfer  for  variations  in  rotation  or  density  ratio  diminishes. 
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7.2.2  Low  Pressure  Surfaces 


The  heat  transfer  from  the  low  pressure  surfaces  from  the  first  and  second 
passages  of  Figures  29  and  30  is  more  complex  than  that  from  the  high  pressure 
surfaces.  Heat  transfer  in  the  first  passage  decreases  with  increasing 
rotation  nvunber  for  low  values  of  rotation  number  (i.e.  fl  d/V  <  0.2  at  the 
downstream  location)  and  then  subsequently  increases  again  with  increases  in 
rotation  for  larger  values  of  rotation  number.  Additionally,  as  with  the  high 
pressure  surfaces  in  the  first  passage,  heat  transfer  increases  with  increases 
in  the  density  ratio.  A  similar  characteristic  in  the  heat  transfer 
distributions  is  observed  in  the  second  passage  for  radially  inflow  as  well. 
However,  with  one  exception,  the  large  effects  of  density  ratio  observed  on 
the  low  pressure  surfaces  of  the  first  passage  are  diminished  in  the  second 
passage.  The  exception  is  that  heat  transfer  is  slightly  increased  with 
increasing  density  ratio  near  the  inlet  of  the  second  passage. 

The  extrapolated  values  of  the  heat  transfer  ratios  as  approaches 

zero  are  modified  from  those  previously  presented  for  the  same  data  set  in 
Wagner  et  al  (1989,  1991)  and  Wagner  et  al  (1990).  Based  on  unpublished  (1990) 
preliminary  results  from  a  mass  transfer  experiment  at  the  University  of 
Darmstadt,  Germany  and  numerical  studies  of  lacovides  and  Launder  (1990)  for 
flows  in  rotating  radial  ducts,  the  authors  now  believe  that  the  heat  transfer 
ratio  of  (  A  P/P) in  =  0  on  the  low  pressure  side  of  the  duct  would  not 
increase  with  increasing  rotation  number  after  reaching  a  minimal  value.  The 
conclusion  from  this  position  is  that  the  heat  transfer  coefficients  for 
conditions  with  (  A  P/P) in  ~  0.2  can  be  as  much  as  3  to  4  times  those 
obtained  when  (  AP/P)in  approaches  zero. 

The  more  complicated  heat  transfer  distributions  on  the  low  pressure  surfaces 
of  the  coolant  passage  are  .attributed  to:  1)  the  combination  of  buoyancy 
forces  and  the  stabilization  of  the  near-wall  flow  for  low  values  of  the 
rotation  number  and  2)  the  developing  Coriolis  driven  secondary  flow  cells  for 
the  larger  values  of  the  rotation  number.  It?  is  postulated  that  the  relatively 
large  effects  from  variations  in  density  ratio  near  the  inlet  of  the  second 
passage  and  the  small  effects  near  the  end  of  the  second  passage  are  due  to 
the  development  of  the  near-wall  thermal  layers.  Near  the  inlet  of  the  second 
passage,  the  thermal  layers  are  postulated  to  be  thin  because  of  the  strong 
secondary  flows  in  the  first  turn  region.  With  increasing  X/d,  the  turn 
dominated  secondary  flows  diminsh  and  the  counteracting  effect  of  buoyancy  and 
the  Coriolis  generated  secondary  flow  increases. 

7.3  Effect  of  Buoyancy  Parameter 

The  analysis  of  the  equations  of  motion  for  flow  in  a  rotating  radial  passage 
by  Suo  (1980)  showed  that  1)  variations  in  the  momentum  of  the  flow  in  the 
plane  perpendicular  to  the  passage  centerline  (cross  stream  flow)  will  be 
proportional  to  the  rotation  number , i) d/V,  and  2)  variations  in  the  moment™ 
of  the  flow  parallel  to  the  passage  centerline  (buoyant  flows)  will  be 
proportional  to  the  buoyancy  parameter  ( A  P/P) (R/d) ( H d/V)2 .  The  buoyancy 
parameter  defined  is  equivalent  to  the  ratio  of  the  Grashof  number  (with  a 
rotational  gravitation  term,  Rfi  to  the  square  of  the  Reynolds  number  and 
has  previously  been  used  to  characterize  the  relative  importance  of  free-  and 
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forced-convection  In  the  analysis  of  stationary  and  mixed-convection  heat 
transfer.  Guidez  (1988)  used  a  similar  analysis  to  establish  appropriate  flow 
parameters  for  the  presentation  of  his  results.  These  parameters,  J2d/V  and 
(AP/P  )(R/d)(nd/V)2,  will  also  be  used  in  the  present  discussion  of  the 
effects  of  Coriolis  and  buoyancy  forces  on  the  heat  transfer. 

The  data  and  extrapolated  results  presented  in  Figures  29  and  30  show  that  the 
effects  of  Coriolis  and  buoyancy  forces  are  coupled  through  the  entire  operat¬ 
ing  range  investigated.  The  results  from  Figure  29  combined  with  those  for  R/d 
=  33,  are  presented  in  Figures  31  and  32  as  the  variation  of  the  heat  transfer 
ratio  with  the  buoyancy  parameter.  The  local  density  ratio  and  radius,  R,  are 
used  in  the  buoyancy  parameter.  For  a  constant  wall  temperature  boundary 
condition,  the  wall-to-bulk  temperature  difference  (hence  the  local  density 
ratio)  decreases  as  the  flow  progresses  downstream,  therefore,  the  range  of 
the  buoyancy  parameter  also  decreases  with  increasing  values  of  X/d.  Results 
for  the  same  value  of  the  rotation  number  are  connected  with  lines  where  the 
results  are  not  well  correlated  by  the  buoyancy  parameter.  The  lines  at 
constant  rotation  number  are  extrapolated  to  the  value  of  the  heat  transfer 
ratio  estimated  for  a  density  ratio  (and  also  buoyancy  parameter)  of  0  as 
described  in  section  7.2. 

The  heat  transfer  ratios  for  this  trailing  side  (Figure  32)  of  the  passage 
increase  with  the  buoyancy  parameter.  The  rate  of  increase  in  the  heat 
transfer  ratio  with  increasing  buoyancy  parameter  is  greatest  at  the  X/d  = 

12.4  location  for  values  of  buoyancy  parameter  less  than  0.4.  For  values  of 
buoyancy  parameter  greater  than  0.4,  the  rate  of  increase  is  less.  Thus,  two 
ranges  of  buoyancy  parameter  appear  to  exist  with  different  heat  transfer 
characteristics.  Generally,  the  heat  transfer  variations  from  the  trailing 
side  form  a  1:1  correspondence  with  the  buoyancy  parameter  (i.e,  singled  value 
function)  and  are  well  correlated  by  the  buoyancy  parameter  for  all  values  of 
X/d  shown. 

Examination  of  the  heat  transfer  results  from  the  leading  side  (Figure  31) 
suggest  that  at  least  three  ranges  of  buoyancy  pareuneter  exist  where  the  heat 
transfer  is  dominated  by  different  fluid  dynamic  mechanisms  (i.e.,  Coriolis, 
buoyancy,  etc.).  At  X/d  =  12.4,  there  is  a  range  of  buoyancy  parameter  for 
values  less  than  0.1  (range  A)  where  the  heat  transfer  ratios  decrease  sharply 
(from  approximately  1.1  to  just  over  0.5)  with  increasing  values  of  the 
buoyancy  parameter.  Within  the  second  range  from  0.1  to  approximately  0.5 
(range  B),  the  heat  transfer  ratios  increase  sharply  with  increasing  values  of 
the  buoyancy  parameter.  For  the  third  range,  with  values  of  the  buoyancy 
parameter  greater  than  0.5  (range  C),  the  heat  transfer  ratio  increases  at  a 
lower  rate,  with  Increasing  values  of  buoyancy  pareuneter.  For  lower  values  of 
X/d,  the  three  ranges  are  less  well  defined.  However,  the  minimum  value  of 
heat  transfer  ratio  (which  defines  the  end  of  the  first  buoyancy  parameter 
range)  occurs  at  increasing  magnitudes  of  buoyancy  parameter  (0.1  to  0.4)  as 
X/d  decreases  from  12.4  to  4.7.  The  heat  transfer  on  the  leading  surface  at 
values  of  X/d  =  4.7  and  8.5  is  governed  by  a  more  complex  relationship  of 
streamwise  distance,  rotation  number  and  buoyancy  parameter.  However,  the 
results  from  the  leading  side  for  X/d  =  12.4  are  well  correlated  by  the 
buoyancy  parameter  for  values  of  the  buoyancy  parameter  greater  than  0.2. 
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Figure  32.-  Variation  of  Heat  Transfer  Ratio  with  Buoyancy  Parameter. 
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The  analyses  of  these  heat  transfer  results  show  that  1)  the  buoyancy 
parameter  correlates  the  heat  transfer  ratio  data  from  the  trailing  side  of 
the  coolant  passage  and  from  the  leading  side  at  the  downstream  location,  2) 
the  data  was  not  correlated  by  the  buoyancy  parameter  near  the  inlet  on  the 
leading  durface  due  to  a  complex  Interaction  of  stabilization  of  the  near  wall 
flow,  buoyancy  forces  and  Coriolis  effects,  and  3)  the  heat  transfer  in 
rotating,  smooth  passages  is  governed  by  complex  interactions  of  the  viscous, 
Coriolis  and  buoyancy  forces  on  the  fluid. 

7.4  Effect  of  Flow  Direction 

The  data  was  also  analyzed  to  determine  the  effects  of  flow  direction 
(radially  inward  or  radially  outward)  on  the  heat  transfer  characteristics  and 
to  determine  the  differences  between  the  first  leg  with  outward  flow 
downstream  of  an  inlet,  the  second  leg  with  inward  flow  downstream  of  an  180“ 
turn  and  the  third  leg  with  outward  flow  downstream  of  an  180“  turn.  The 
variations  of  the  heat  transfer  ratio  with  buoyancy  parameter  for  the  heater 
test  section  most  downstream  of  the  inlet  or  turn  for  each  of  the  three  legs 
are  presented  in  Figures  33  and  34.  This  is  the  streamwise  location  in  each 
leg  at  which  the  asymptotic  heat  transfer  characteristics  were  approached  for 
lid/V  =  0,  (see  Figure  15). 

The  data  presented  in  Figures  29  and  30  showed  that  the  effects  of  Coriolis 
and  buoyancy  forces  are  coupled  in  the  first  two  passages  through  the  entire 
operating  range  investigated.  The  results  from  Figures  29  and  30  plus  addition¬ 
al  results  from  the  third  passage  are  combined  with  those  for  R/d  =  33  and  are 
presented  in  Figures  33  and  34  as  the  variation  of  the  heat  transfer  ratio 
with  the  buoyancy  parameter  based  on  the  local  density  ratio  and  radius,  R. 

Note  that  in  Figures  33  and  34,  the  range  of  the  buoyancy  parameter  data 
becomes  more  compressed  as  the  fluid  progresses  downstream.  This  is  because 
the  wall-to-bulk  temperature  difference  (hence  the  local  density  ratio)  is 
decreasing  as  the  fluid  progresses  downstream  and  the  bulk  temperature 
increases.  The  ranges  of  heat  transfer  ratio  for  the  last  location  in  the 
first  passage  is  shown  as  a  shaded  band  with  the  results  from  the  second  and 
third  passages  for  comparison. 

Heat  transfer  distributions  from  the  low  pressure  surfaces  of  each  of  the 
three  passages  (Figure  33)  exhibit  a  decrease  with  increasing  values  of 
buoyancy  between  0.0  andO.15.  Heat  transfer  subsequently  increases  again  with 
increasing  values  of  buoyancy.  Heat  transfer  on  the  low  pressure  surfaces  of 
rotating  coolant  passages  is  governed  by  complex  relationships  of  streamwise 
location,  rotation  number,  and  buoyancy  parameter.  However,  the  heat  transfer 
results  are  reasonably  well  correlated  in  the  first  two  passages  by  the 
buoyancy  parameter  for  values  of  buoyancy  parameter  greater  than  0.2. 

The  heat  transfer  results  from  the  high  pressure  surfaces  in  the  first  passage 
(Figure  34)  are  correlated  well  by  the  buoyancy  parameter.  The  second  passage 
with  radially  inward  flow  had  different  heat  transfer  characterisitcs  than  the 
first  and  third  passages  with  radially  outward  flow.  Whereas  the  heat  transfer 
ratios  for  the  high  pressure  surfaces  of  the  first  and  third  passages 
increased  with  the  buoyancy  parameter,  the  heat  transfer  in  the  second  passage 
was  lower  and  relatively  independent  of  buoyancy  parameter  for  values  of 
buoyancy  greater  than  0.05.  These  results  for  co-flowing  and  counter-flowing 
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Figure  33.-  Comparison  of  Heat  Transfer  Ratios  from  Passages  with  Inward  and 
Outward  Flow. 
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buoyancy  effects  on  the  high  pressure  surfaces  are  generally  consistent  with 
the  stationary  combined  free-and-forced  convection  experiments  of  Eckert  et 
al.  (1953).  They  measured  decreased  levels  of  heat  transfer  for  the  co-flowing 
conditionC  i.e.  similar  to  that  of  radially  inward  flow  in  rotating  sytems). 

7.5  Conclusions 

The  analysis  of  the  experimental  results  to  determine  the  combined  effects  of 
forced  convection,  buoyancy  and  Coriolis  forces  on  heat  transfer  in  smooth 
rotating  radial  passages  has  produced  several  interesting  and  remarkable 
conclusions : 

1.  Both  the  density  ratio, (  Z\P/P)iri»  the  rotation  number,  Qd/V  were 
found  to  cause  large  changes  in  the  heat  transfer  ratio,  as  much  as 
factors  of  2.5  and  3.5,  respectively. 

2.  The  flow  field  was  never  completely  dominated  by  free  convection  for  the 
test  conditions  in  this  study.  The  conclusion  is  easily  discerned  by 
observing  the  variation  of  the  heat  transfer  ratio  between  the  leading 
and  trailing  surfaces. 

3.  The  heat  transfer  ratio  was  found  to  be  primarily  a  function  of  the 
buoyancy  parameter  on  the  low  pressure  surfaces  (i.e.,  leading  surfaces 
for  outward  flow  and  trailing  surfaces  for  inward  flow)  for  values  of 
the  buoyancy  parameter  greater  than  0.2  and  for  X/d  =  12.4. 

4.  The  heat  transfer  ratio  was  found  to  be  primarily  a  function  of  the 
buoyancy  parameter  for  the  high  pressure  side  of  the  passage  for  flow 
outward  only  (i.e.,  trailing  side  of  passage). 

5.  The  variations  of  heat  transfer  ratio  with  buoyancy  parameter  were 
approximately  the  same  for  the  furthest  downstream  straight  section 
segment  on  the  low  pressure  side  of  all  three  passages.  Buoyancy  had 
essentially  the  same  effect  on  heat  transfer  on  the  low  pressure  side 
whether  the  flow  is  inward  or  outward. 

6.  The  heat  transfer  ratio  on  the  high  pressure  side  of  the  passage  is 
significantly  affected  by  flow  direction.  The  effects  of  buoyancy  and 
Coriolis  forces:  (a)  combined  to  increase  heat  transfer  on  the  flow 
outward  trailing  surfaces;  and  (b)  cancelled  on  the  flow  inward  leading 
surfaces  with  the  result  that  the  heat  transfer  ratio  is  approximately 
constant.  Possible  explanations  for  this  phenomena  include  changes  in 
the  turbulence  structure  due  to  combinations  of  flow  direction  and 
buoyancy  force  direction. 

7.  The  heat  transfer  ratios  for  R/d  =  49  and  33  are  well  correlated  by  the 
buoyancy  parameter  at  the  downstream  test  surface  of  each  passage. 
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8.0  CORRELATION  OF  HEAT  TRANSFER  RESULTS 


8.1  Physical  Effects  of  Rotation 

It  is  well  docvimented  that  secondary  flow  circulations  and  shear  layer 
stability  variations  are  the  principal  manifestations  of  rotated  channel 
flows.  These  effects  are  illustrated  in  Figure  35  which  depicts  a  rotated 
square  duct  with  radially  outward  flowing  coolant.  The  Coriolis  force 
accelerates  the  low  momentum  fluid  on  the  sidewalls  toward  the  low  pressure 
(leading)  surface  where  the  boundary  layer  is  stabilized,  in  some  cases, 
becoming  laminar.  The  core  fluid  shifts  toward  the  high  pressure  (trailing) 
surface  where  boundary  layer  is  highly  destabilized,  undergoing  a  turbulent 
bursting  process  similar  to  that  found  on  concave  surfaces  (Taylor-Goertler 
roll  cells).  It  is  important  to  understand  that  the  flow  on  each  of  the  three 
surface  types  (leading,  trailing,  and  sidewalls)  in  the  square  duct 
experiences  a  different  set  of  physical  changes  during  rotation. 

Moore  (1967)  measured  the  centerline  velocity  profile  in  a  long  rotating 
square  duct  shown  in  Figure  36.  Moore's  data,  shown  in  Figure  37,  validates 
the  shifting  of  the  core  fluid  from  the  center  of  the  duct  to  the  high 
pressure  surface.  The  boundary  layer  becomes  thick  on  the  low  pressure  surface 
and  thin  on  the  high  pressure  surface.  In  a  heated  duct,  this  effect  would 
insulate  the,  low  pressure  wall  from  the  cool  core  fluid  and  move  this  core 
fluid  closer  to  the  high  pressure  wall,  thus  Increasing  the  heat  pickup  there. 

Moon  (1964)  examined  the  affect  of  rotation  on  internal  flow  in  a  rectangular 
duct  shown  in  Figure  38.  He  measured  centerline  velocities  at  five  locations 
in  183  cm  (72")  long  rotating  duct.  These  are  plotted  along  with  the  boundary 
layer  and  displacement  thickness  in  Figure  39.  The  velocity  profiles  and 
boundary  layer  thickness  both  show  that  the  low  pressure  side  (suction) 
boundary  layer  becomes  very  thick  as  the  flow  moves  downstream.  At  the  same 
time,  the  boundary  layer  on  the  high  pressure  side  of  the  duct  maintains  the 
classical  turbulent  boundary  layer  shape,  except  at  the  end  of  the  duct,  where 
it  becomes  thin.  This  shifting  of  the  boundary  layers  on  both  surfaces  occurs 
as  the  Coriolis  acceleration  transports  fluid  from  the  sidewalls  onto  the  low 
pressure  surface.  Note  that  this  thickening  on  the  low  pressure  side  is  a 
strong  function  of  distance  from  the  duct  inlet. 

Using  previously  measured  skin  friction  coefficients  and  the  Reynolds  analogy, 
we  can  generally  predict  the  effects  of  the  Coriolis  force  on  local  heat 
transfer.  Moore  (1967)  shows  in  Figure  40  that  at  two  rotation  ntunber 

conditions,  the  skin  friction  on  the  high  and  low  pressure  surfaces  varies 

significantly  and  is  a  strong  function  of  channel  aspect  ratio.  Moore's  data 

at  unity  aspect  ratio  (square  duct)  shows  that  the  low  pressure  surface  skin 

friction  decreases  to  65  percent  of  the  stationary  value.  A  straight  Reynolds 
analogy  would  predict  local  heat  to  decrease  by  35  percent  at  these 
conditions.  This  magnitude  of  decrease  is  realized  in  the  present  square  duct 
experiment,  but  the  decrease  occurs  for  higher  rotation  number  conditions; 
most  likely  since  the  duct  in  this  experiment  is  relatively  short  and  has  less 
length  available  for  Coriolis  secondary  flow  development. 
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Figure  39.-  Moon's  Centerline  Velocity  Profiles  with  Boundary  Layer  and 
Displacement  Thicknesses  (Moon,  1964). 
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Figure  40.-  Moore’s  Skin  Friction  Variation  with  Aspect  Ratio  (Moore,  1967). 


The  aforementioned  results  of  Moore  and  Moon  exeunine  the  effects  of  rotation 
on  non-heated  channel  flows.  The  principle  effect  is  the  generation  of 
secondary  flow  by  the  Coriolis  acceleration  and  the  resulting  wall  shear  layer 
modifications.  In  the  case  of  rotating  heated  flows,  centrifugal  buoyancy 
forces  become  an  important  factor. 

Figure  41  deplicts  the  directions  of  the  convective  and  buoyancy  forces  during 
rotation.  Certrifugal  buoyancy  forces  accelerate  the  cold,  more  dense  fluid 
particles  away  from  the  center  of  rotation,  while  the  hot  and  lighter  fluid 
particles  tend  toward  the  axis  of  rotation.  In  this  experiment  the  hot  fluid 
along  the  heated  walls  tends  toward  the  center  of  rotation;  in  the  first  and 
third  passages  of  the  model,  this  direction  opposes  that  of  the  mainstream 
flow.  In  the  second  model  passage,  where  the  coolant  flows  toward  the  axis  of 
rotation,  the  mainstream  and  buoyancy  force  directions  become  aligned. 

Eckert  (1954)  defined  these  two  buoyancy  conditions  as  "counterflow"  and 
"parallel  flow"  and  examined  them  in  detail  in  nonrotating  experiments.  In 
Eckert’s  case,  gravity  rather  than  certrifugal  forces  provides  the  buoyancy 
acceleration.  In  a  later  section,  his  results  will  be  compared  to  the  data  of 
this  experiment.  Eckert’s  results  help  explain  some  of  the  different  trends 
seen  in  radially  inward  versus  radially  outward  flowing  coolant  passages. 

Buoyancy,  therefore,  complicates  the  three  dimensional  flow  established  by 
Coriolis  secondary  flow.  The  heated  wall  shear  layer  always  tends  toward  the 
axis  of  rotation  and  Interacts  with  the  cross-stream  secondary  flows  generated 
by  the  Coriolis  force.  In  the  smooth  square  duct,  the  result  is  a  heat 
transfer  distribution  strongly  dependant  on  local  conditions,  both  axially  and 
circumferentially.  Boundary  layers  on  each  of  the  three  surface  types 
(leading,  trailing  and  sidewalls)  are  influenced  differently  by  the  forces  of 
rotation. 


66 


RADIAL 

OUTFLOW 

PASSAGE 

length 


FLOW 

INLET 


Figure  41.-  Mainstream  Flow  Direction  and  Buoyancy  Force  Direction  in  Radially 
Inward  and  Radially  Outward  Flowing  Passages. 

8.2  Rotation  Heat  Transfer  Correlations 

In  this  section,  each  passage  region  (sidewall,  high  pressure  wall,  turn,  and 
low  pressure  wall)  will  be  correlated  in  a  manner  that  is  consistent  with  the 
dominate  physics  for  that  particular  region. 

8.2.1  Side  Wall  Heat  Transfer 

This  section  covers  the  correlations  of  the  sidewall  heat  transfer  data  in  the 
first  radially  outward  flowing  passage.  These  correlations  are  presented  first 
because  the  side  wall  heat  transfer  is  the  best  behaved  during  rotation.  On 
these  surfaces,  the  heat  transfer  is  dominated  by  the  development  of 
Coriolis-induced  circulations  and  influenced  less  than  the  leading  and 
trailing  surfaces  by  the  stabilizing/destabilizing  effects  of  rotation. 

The  side  wall  heat  transfer  results  for  Re  =  25,000  for  the  first  straight 
passage  are  shown  in  Figure  42.  These  data  correspond  with  test  section 
surfaces  numbered  B,  C  and  D  shown  in  Figure  8.  The  data  is  plotted  as  the 
ratio  of  rotating  to  stationary  (measured)  Nusselt  number  versus  the 
rotational  Grasshof  number  divided  by  a  Reynolds  number  squared.  This  ratio 
breaks  down  into  the  following  combination  of  the  basic  dimensionless 
parameters ; 


8.1 
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LEG  1  -  FLOW  OUTWARD  Re  “  25,000,  R/d  =  49,  a  =  0 
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Figure  42.-  Sidewall  Heat  Transfer  at  Re  =  25,000  for  the  First  Straight 
Passage . 

Each  of  these  basic  ratios  effects  the  rate  of  heat  transfer  in  rotating 
internal  flows.  The  first  parameter  on  the  left-hand  side  of  the  equation  is 
the  rotation  number  Ro  =nd/V;  it  equals  the  inverse  of  the  Rossby  number.  The 
rotation  n\imber  is  a  ratio  of  the  Coriolis  force  to  the  inertia  force.  The 
second  term  represents  the  nondimens ional  radius  ratio.  The  third  term  is  an 
entry  length  or  development  length  ratio  and  the  last  is  a  temperature  or 
density  ratio. 

This  method  of  plotting  the  data  reveals  an  important  trend.  The  heat  transfer 
behavior  changes  with  increasing  rotation  number.  To  facilitate  correlating 
these  trends,  the  data  is  separate  into  two  groups:  Regimes  of  low  and  high 
rotation  ntimber.  In  Figure  42,  the  steepness  in  slope  of  the  curves  of 
constant  rotation  number  steadily  increases  with  increases  in  the  buoyancy 
par^eter;  Grashoff  number  divided  by  Reynolds  number  squared.  At  low  rotation 
numbers,  the  curves  remain  relatively  flat  while  for  the  higher  rotation 
numbers  the  slopes  increase.  Below  and  above  Ro=0.20,  the  rotation  and 
buoyancy  forces  interact  differently  and  correlate  with  different  parameters. 
It  should  be  noted  that  this  rotation  number  level  Ro  =  0.20,  will  also  be 

shown  to  mark  the  beginning  of  important  trends  on  other  surfaces  in  the 
passage . 
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a)  Side  Wall  Heat  Transfer  Low  Rotation  Regime  —  Figure  43  shows  the 
sidewall  heat  transfer  correlated  for  the  low  rotation  regime  Ro  <  0.20.  The 
correlating  parameter  is  a  product  of  the  rotation  number  and  the  buoyancy 
parameter  defined  as  the  rotational  Grasshof  number  divided  by  the  Reynolds 
number  squared.  The  resulting  parameter  breaks  down  as  follows; 


The  low  rotation  regime  sidewall  heat  transfer  data  may  be  correlated  as 
follows: 


assuming  Rq  0-5(R/d)(X/d)C\T/T)<0.9;  Nu/NUq  =1.0 
and  for  Rq  0’5(8/d)(X/d)(AT/T)>0.9;  the  following  equation  holds 


Nu 

Nu^ 
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0.6351  Ro  “  5 


ft  X 

d  d  T  I 


(8.3) 


This  equation  simplifies  into  the  form 


Nu 

Wn 


=  0.635  Ro  °  (  - 


(;  s  f  r 


(8.4) 


b)  Side  Wall  Heat  Transfer;  Hieh  Rotation  Regime  —  Using 
(Grx/Re2x)Ro-0-8  =  Rol-2(R/d)(X/d)(  At/T)  as  the  correlation  parameter, 
Figure  44  presents  the  side  wall  heat  transfer  correlation  for  the  high 
rotation  regime,  Ro  2.  0.20.  The  data  collapses  reasonably  well  for  the 
combination  of  parameters  shown,  all  but  a  few  points  fall  within  10%  of  the 
equation; 


(8.5) 

The  equation  reveals  the  relative  strength  of  each  of  the  individual 
dimensionless  parameters  on  rotating  heat  transfer  for  the  side  wall. 

c)  Side  Wall  Heat  Transfer;  Reynolds  Number  Effect —  In  Figures  43  and  44, 
the  effects  of  Reynolds  number  variations  are  also  evaluated.  The  tests 
conducted  where  the  Reynolds  number  varied  from  the  base  condition  of  Re  = 
25,000  are  indicated  on  the  figures  as  half  open  symbols.  The  open  symbols 
represent  Re  =  25,000. 

In  Figure  43  the  test  conducted  at  Re  =  50,000  and  shown  as  the  only  half  open 
symbols  on  the  figure  agree  very  well  with  this  method  of  correlating  the  low 
rotation  number  regime  data. 
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LEG  1  -  FLOW  OUTWARD,  a  =  0 
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Figure  44.-  Effect  of  Radius  Variation  on  Sidewall  Heat  Transfer  for  High 
Rotation.  (Ro  >  0.20). 


Figure  44,  correlating  the  high  rotation  regime  dat\,  compares  both  high  and 
low  Reynolds  number  cases  with  the  data  for  Re  =  25,000.  Test  number  5  at 
Re  =  12,600  and  test  number  6  at  50,000  were  both  conducted  at  the  same 
rotation  number,  Ro  =  0.25.  The  higher  Reynolds  number  data  agrees  very  well 
with  the  previously  correlated  data  except  for  a  single  point.  The  lower 
Reynolds  number  data  appears  slightly  higher  than  the  correlation  curve,  yet 
still  remains  less  than  15%  above  this  curve.  It  is  believed  that  low  Reynolds 
number  flow  is  more  severely  affected  by  rotation.  For  the  moderate  to  high 
Reynolds  number  flows  tested,  the  buoyancy  parameter  adequately  accounts  for 
the  effects  of  Reynolds  number  variations. 
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d)  Side  Wall  Heat  Transfer  Model  Radius  Effect  —  All  the  open  and  half  open 
symbol  data  correlated  to  this  point  has  been  for  the  nondimensional  model 
radius  ratio  R/d  =  49.  Four  tests  were  conducted  at  a  smaller  radius,  R/d  = 

33.  Figure  44  presents  this  data  as  solid  symbols.  The  behavior  is  similar  to 
that  for  the  larger  model  radius.  The  buoyancy  parameter,  therefore,  accounts 
for  radius  ratio  effects  on  heat  transfer  for  the  range  tested. 

8.2.2  High  Pressure  Wall  Heat  Transfer 

This  section  correlates  the  heat  transfer  results  for  the  first  passage  on  the 
high  pressure  side  of  the  smooth  square  channel.  This  is  a  radially  outward 
flowing  passage,  and  the  high  pressure  surface  wall  sections  correspond  to 
those  numbered  50,  51,  and  52  on  the  trailing  side  of  the  model  shown  in 
Figure  8. 

The  heat  transfer  results  are  correlated  using  a  combination  of  the  basic 
nondimensional  parameters  shown  in  the  previous  section. 


9^ 

V 


(8.6) 


In  Figure  45  the  heat  transfer  data  for  Re  =  25,000  and  radius  ratio  R/d  =  49 
are  shown.  The  data  are  plotted  as  the  ratio  of  rotating  to  stationary  heat 
transfer  versus  the  nondimensional  rotational  buoyancy  parameter,  Gr^/Re^j^. 
This  parameter  breaks  down  into  the  group  of  dimensionless  parameters  listed 
above,  where  m=2  and  the  powers  n,  p  and  q  are  all  equal  to  one.  Lines  of 
constant  rotation  number,  Ro,  connect  wall  test  sections  numbered  50,  51,  and 
52  in  Figure  8. 

Plotted  in  this  manner,  the  data  reveal  an  Important  trend.  As  the  rotational 
buoyancy  parameter  increases,  and  the  rotation  number  surpasses  the  critical 
level  Ro  =  0.20,  the  data  tend  toward  falling  on  a  uniform  curve.  But  at  lower 
levels  of  this  parameter,  there  exist  distinct  curves  of  constant  rotation 
nmber.  It  is  believed  that  centripetal  buoyancy  forces  become  more  dominant 
at  higher  rotation  rates.  As  would  be  expected,  the  temperature  or  density 
ratio  directly  affects  the  behavior  of  the  heat  transfer  during  rotation.  The 
conclusion  reached  here  is  that  the  complicated  interaction  between  Coriolis 
and  centripetal  buoyancy  forces  greatly  changes  for  low  and  high  rotation 
number  ranges.  As  a  result,  the  following  figures  correlate  the  data  according 
to  the  two  distinct  groups  mentioned  previously;  low  rotation  and  high 
rotation  number  regimes. 
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LEG  1-  FLOW  OUTWARD  Re^25.000  ,  R/d  =  49  ,  a  =  0 
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Figure  45.- 


High  Pressure  Surface  Heat  Transfer  for  Radially  Outward  Flow 
Re  =  25,000. 


a)  High  Pressure  Side:  Low  Rotation  Number  Regime  Rp  <  0.20  —  Heat 
transfer  on  the  high  pressure  side  of  the  outward  flowing  passage  for  the 
rotation  number  regime  is  correlated  on  Figure  46.  The  data  correlates 
reasonably  well  with  the  parameter 
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The  open  symbols  represent  data  where  radius  and  Reynolds  number  are  held 
constant  and  rotation  rates  and  temperature  ratios  are  varied.  Also  X/d  is 
varied  as  data  is  acquired  for  3  locations  (Sections  50,51,52  of  Fig.  8)  in 
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Figure  46,-  Effect  of  Reynolds  Number  for  Low  Rotation  (Ro  <  0.20)  on  High 
Pressure  Surface  Heat  Transfer. 


the  straight  section.  The  half  open  symbols  indicate  two  higher  Reynolds 
number  test  conditions.  The  largest  discrepancies  occur  for  those  test 
sections  located  farthest  downstream  of  the  entrance  or  just  before  the  turn. 
Here  the  correlation  overpredicts  the  data.  It  is  suspected  that  the  turn 
effects  propagate  upstream  and  affect  these  few  data  points.  In  general  the 
low  rotation  number  regime  data  are  well  correlated  by  the  above  equation. 

b)  High  Pressure  Side:  High  Rotation  Number  Regime  Ro>0.20  —  Figure  47 
correlates  the  high  pressure  surface  heat  transfer  data  for  rotation  number  Ro 
^0.20.  This  correlation  incorporates  singular  variations  in  rotation  rate, 
Reynolds  number,  radius,  density  or  temperature  ratio,  and  location  in  the 
passage.  The  data  best  collapses  with  the  parameter 
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LEG  1  -  FLOW  OUTWARD,  a  =  0 


8.2.3  Turn  Heat  Transfer 


The  experimental  model  shown  in  Figure  8  includes  three  instrxunented  turns. 

The  first  and  third  turns,  referred  to  as  the  tip  turns,  behave  very 
differently  from  the  root  turn  (second  turn)  during  rotation.  Centripetal 
buoyancy  drives  the  rotating  turn  heat  transfer  and  the  parameter 
(iiR/V)(AT/T)  correlates  the  data.  Each  turn  contains  a  total  of  eight  heated 
wall  sections  where  heat  transfer  coefficients  were  obtained.  The  inside  wall 
of  each  turn  was  not  Instrumented. 

a)  Tip  Turn  Heat  Transfer  —  Figure  48  presents  the  tip  turn  outside  wall 
heat  transfer  results  for  all  values  of  Ro  and  AT/T  tested  at  Re  =  25,000. 
Plotted  versus  the  centripetal  buoyancy  parameter,  ((iR/V) (AT/T) ,  the  heat 
transfer  ratio  is  correlated  at  each  wall  section  or  streamwise  location.  The 
largest  increases  occur  at  the  beginning  of  the  turn  and  gradually  decrease 
through  the  turn.  Solid  symbols  correspond  to  wall  sections  5  and  7,  and  open 
symbols  correspond  to  sections  6  and  8. 

Figure  49(a  and  b)  show  the  leading  (open  symbols)  and  trailing  (solid 
symbols)  surface  heat  transfer  results  plotted  with  the  centripetal  buoyancy 
parameter.  The  upper  portion  of  the  figure  examines  the  first  half  of  the  turn 
and  the  lower  portion  shows  the  second  half.  In  the  first  half  of  the  turn, 
the  trailing  wall  heat  transfer  is  the  higher  of  the  two.  As  the  flow  turns 
the  corner,  reversing  from  radially  outward  to  radially  inward,  the  leading 
surface  heat  transfer  becomes  the  higher  of  the  two  sides  in  the  second  half 
of  the  turn. 

Unlike  the  results  of  the  straight  sections,  the  heat  transfer  on  both  the 
leading  and  trailing  surfaces  of  the  tip  turn  is  generally  augmented  with 
rotation.  The  low  pressure  surface  does  not  experience  a  decrease  in  heat 
transfer  because  the  centripetal  buoyancy  forces  dominate  the  Coriolis  forces. 
As  the  flow  negotiates  the  turn,  the  cross  channel  pressure  gradient  drives 
the  low  momentum  fluid  toward  the  inside  of  the  turn.  Rotation  causes  the  less 
dense  wall  fluid  to  gravitate  toward  the  inside  of  the  turn  as  the  denser 
mainstream  fluid  accelerates  away  from  the  axis  of  rotation.  In  the  tip  turn 
the  centripetal  buoyancy  forces  enhance  the  secondary  motions  already  generated 
by  the  curvature  of  the  turn;  the  secondary  flows  resulting  from  both  curvature 
and  rotation  are  aligned  parallel  and  complement  one  another.  The  result  in 
the  tip  turn  is  an  intensification  in  the  secondary  motions  normally  generated 
in  a  turn  and  significant  increases  in  the  heat  transfer  during  rotation. 

b)  Root  Turn  Heat  Transfer  —  All  surfaces  in  the  root  turn  experience  a  drop 
in  heat  transfer  during  rotation.  Similar  to  the  tip  turn,  the  centripetal 
forces  accelerate  the  dense  mainstream  fluid  away  from  the  axis  of  rotation. 

In  the  root  turn,  this  direction  is  toward  the  inside  wall  of  the  turn.  The 
mainstrean  coolant  that  normally  impinges  on  the  outside  wall  as  the  flow 
turns  is  acted  upon  by  rotation  and  inhibited  from  generating  stream  secondary 
flows.  As  a  result,  the  heat  transfer  enhancing  mechanism  of  the  root  turn  is 
hindered. 
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Figure  49.-  Leading  and  Trailing  Surface  Heat  Transfer  in  Tip  Turn 
(Solid  Symbols  =  Trailing  Surface). 


It  is  interesting  to  note  that  during  rotation  the  root  turn  heat  transfer 
almost  reduces  to  the  level  of  heat  transfer  in  a  straight  duct.  This  suggests 
that  the  secondary  flows  in  the  root  turn  can  almost  be  eliminated  by 
centripetal  forces.  For  example,  at  wall  section  no.  27  (Figure  50),  the  heat 
transfer  (Nu/NUq=  0.41)  decreases  from  the  stationary  NUq  =  165  to  the 
rotating  Nu  =  68  at  rotation  number,  Ro  =  0.24  (test  no.  4).  This  drop  in  heat 
transfer  goes  from  315%  to  130%  of  stationary  fully  developed  turbulent  heat 
transfer  (Nu^  =  52.3).  Essentially,  the  centripetal  buoyancy  forces  reduce 
the  root  turn  secondary  flows  to  the  magnitude  where  the  flow  turns  the  comer 
as  if  it  was  in  a  straight  duct. 

Figure  50  illustrates  the  root  turn  results  for  the  outside  walls.  As  noted, 
rotation  decreases  heat  transfer  everywhere  in  the  root  turn.  The  upper 
portion  of  Figure  50  presents  the  wall  sections  numbered  26  and  29  which  are 
located  at  the  beginning  and  end  of  the  turn,  while  those  numbered  27  and  28, 
located  in  the  bucket  region  of  the  turn,  are  shown  in  the  lower  portion.  They 
are  presented  in  this  manner  because  of  their  similar  behavior  with  rotation. 
The  two  sections  in  the  bucket  region  shown  the  larger  decrease  in  heat 
transfer,  leveling  off  near  45%.  The  beginning  and  end  of  the  turn  drop  to  60% 
of  the  stationary  levels. 

Figure  51  quantifies  the  decreases  in  root  turn  heat  transfer  on  the  leading 
and  trailing  surfaces.  All  four  of  these  surfaces  settle  around  60%  of  the 
stationary  measured  heat  transfer  rates. 

c)  Conclusion;  Turns  —  The  rotation  of  internally  cooled  turns  significantly 
changes  the  heat  transfer  from  the  non~rotated  levels.  Generally,  non-rotated 
turns  show  a  200-300%  increase  over  fully  developed  straight  duct,  turbulent 
heat  transfer.  In  rotated  tip  turns,  these  levels  are  further  increased  and  in 
root  turns,  they  are  reduced.  The  dominating  mechanism  is  centripetal  buoyancy 
which  sets  up  forces  either  opposed  to,  or  in  the  same  direction  as,  the  turn 
induced  secondary  flows. 

8.2.4  Low  Pressure  Wall  Heat  Transfer 

This  section  examines  the  heat  transfer  on  the  leading  surface  in  the  first 
radially  outward  flowing  passage  of  the  model.  During  rotation  this  surface 
becomes  the  low  pressure  surface  as  Coriolis  forces  accelerate  low  momentum 
fluid  on  the  sidewalls  toward  this  surface.  The  isolated  effect  of  rotation 
number  for  Reynolds  number  of  Re  =  25,000  is  depicted  in  Figure  52.  As  the 
rotation  number  increases,  the  heat  transfer  ratio  decreases  along  the  passage 
near  the  inlet.  For  all  of  the  remaining  locations  on  the  leading  side  of  the 
passage,  the  heat  transfer  ratio  decreases  and  then  increases  again  with 
increasing  rotation  number.  Examination  of  the  leading  side  results  shows  that 
the  location  of  the  local  minimum  in  the  heat  transfer  ratio  for  each  rotation 
number  moves  toward  the  inlet  of  the  passage  as  the  rotation  number  is 
Increased. 

Figure  53  plots  all  the  leading  side  data  for  At  =  44.4°C  (80°F)  and 
Re  =  25,000.  The  data  is  plotted  as  Nu^,  X  being  the  distance  from  the 
inlet,  versus  a  rotational  Rayleigh  number.  The  Rayleigh  number  was  chosen  as 
a  correlation  parameter  in  order  to  determine  what  flow  regime  (laminar, 
transitional  or  turbulent)  the  low  pressure  wall  data  was  in. 
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Figure  51.-  Leading  and  Trailing  Surface  Heat  Transfer  in  Root  Turn 
(Solid  Symbols  =  Trailing  Surface). 
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Figure  53.-  Leading  Side  Heat  Transfer  for  First  Passage:  Radially  Outward 
Flow  (Test  section  elements  34,  35,  36). 


Note  for  each  line  of  constant  Ro  there  are  three  data  points.  Each  of  these 
corresponds  to  one  of  the  three  test  section  elements  downstream  of  the  guard 
heaters  at  the  inlet:  elements  34,  35,  and  36  in  Figure  8. 

For  the  higher  rotation  rates,  Ro  x  0.25,  the  heat  transfer,  plotted  as  Nu^, 
tends  to  collapse  on  a  single  curve.  This  curve  attains  a  minimum  around 
Nujj  =  210  and  begins  a  sharp  upturn  at  Ra^  =  10^^.  Based  on  the  results 
to  date,  it  is  believed  that  for  the  lower  rotation  rates,  the  data  are 
predominantly  governed  by  Coriolis  forces,  while  at  the  highest  rates 
centripetal  buoyancy  dominates.  There  is  a  flow  regime  between  these  two 
extremes  where  the  wall  shear  layer  is  believed  to  be  laminar. 

Figure  54  plots  Nu^  versus  Ra^,  for  high  rotation,  Ro  x  0.25,  and  it 
includes  all  four  temperature  cases  (Figure  53  included  only  AT  =  44.4®C 
(80°F).  Two  important  points  should  be  emphasized.  Firstly  for  10^®  <  Ra^ < 
10^^  the  heat  transfer  is  constant  at  Nujj  =  210.  This  level  defines  the 
minimum  heat  transfer  attained  for  Re  =  25,000.  Secondly,  for  Ra^  >  10^^, 

Nux  increases  significantly.  This  increase  for  large  Ra^  is  believed  to  be 
induced  by  centripetal  buoyancy  forces. 
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Figure  54.-  Leading  Side  Heat  Transfer  Ro  i_  0.25:  Re  =  25,000 
(Test  section  elements  34,  35,  36). 


Consider  the  range  of  data  where  the  heat  transfer  is  constant  at  Nu^  =  210 
(Figure  54).  If  this  data  is  compared  to  both  the  fully  turbulent  stationary 
heat  transfer  and  to  Kays  (1966),  analytical  solution  for  laminar  flat  plate 
heat  transfer:  Nujj  =  0.565  Pr^^^  Re^^^/^  (Figure  55),  it  appears  to  be 
nearly  identical  to  the  laminar  correlation,  thus  supporting  the  hypothesis 
that  a  flow  regime  containing  laminar  shear  layer  does  exist. 

These  results  can  be  further  substantiated  by  examining  the  work  of  J.P. 
Johnston  (1970).  In  his  rotating  channel  experiment  at  Stanford  University,  he 
discovered  regions  where  the  boundary  layer  on  the  leading  wall  was  completely 
laminar  for  Reynolds  numbers  as  high  as  Re  =  15,000.  Figure  56  schematically 
depicts  the  important  characteristics  of  this  flowfield. 

On  the  leading  side  of  the  7:1  aspect  ratio  channel,  Johnston  saw  a  nearly 
parabolic  mean  velocity  profile  and  an  absence  of  the  bursting  process 
normally  seen  at  the  wall  in  turbulent  flow.  On  the  trailing  side  of  the 
channel  Taylor-Goertler  type  vortices  developed  within  the  turbulent  section 
of  the  mean  velocity  profile.  The  laminar  boundary  layer  on  the  leading 
surface  and  the  highly  turbulent  boundary  layer  with  longitudinal  vortices  on 
the  trailing  surface  help  explain  the  large  decreases  and  increases  seen  in 
heat  transfer  with  rotation. 
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One  of  Johnston's  conclusions  was  that  the  rotation  induced  re-laminarization 
was  highly  Reynolds  number  dependent.  To  evaluate  this  dependency,  the  present 
data  was  examined  at  varying  Reynolds  numbers. 

Figures  57  and  58  compare  the  heat  transfer  results  for  both  rotating 
(Roa;0.24)  and  non-rotating  test  conditions  at  Re  =  12,500  and  Re  =  50,000, 
respectively,  to  Kays  correlation  for  laminar  flow.  The  lower  Reynolds  number 
flow  case.  Figure  57,  matches  the  Kays  correlation  (Kays,  1966)  at  test 
sections  34  and  35,  while  the  higher  Reynolds  number  flow  case  (Figure  58) 
never  reaches  the  laminar  level,  although  it  approaches  this  minimum.  This 
data  indicates  that  at  high  rotation  rates,  the  boundary  layer  on  the  leading 
wall  is  more  likely  to  be  laminar  at  lower  Reynolds  niunber.  This  work  compares 
very  well  with  Johnston's  results. 

Johnston  used  flow  visualization  techniques  to  establish  when  the  leading  wall 
would  laminarize.  Figure  59  extends  Johnston's  mapping  of  flow  regimes  with 
rotation.  Data  at  three  different  Reynolds  numbers  (12,600,  25,200  and  49,600 
for  Ro=0.24)  is  shown  to  complement  Johnston's  results.  At  the  low  Reynolds 
niomber  the  leading  side  is  laminar.  At  the  high  Reynolds  number  the  shear 
layer  is  most  likely  transitional;  somewhere  close  to  but  not  yet  laminar.  At 
the  middle  Reynolds  number  flow.  Re  =  25,200,  there  exists  extensive  data  to 
clarify  where  the  flow  becomes  laminar.  Remember  this  region  existed  when 
NUjj  =  210  (Figure  54).  This  is  when  10^®  <  Ra^  <  lO^^,  indicating  that 
buoyancy  plays  an  important  role  in  the  laminarization  process.  Two  important 
results  need  to  be  highlighted: 

1.  For  the  cases  examined,  the  minimum  level  of  heat  transfer  attained  is 
predicted  by  Kays'  laminar  flat  plate  correlation. 

2.  Centripetal  buoyancy  limits  this  laminarization  process.  For 

Re  =  25,000  and  Ra^  >  10^^  the  leading  side  shear  layer  becomes 
transitional  and  the  heat  transfer  increases. 

Thus  the  low  pressure  surface  wall  shear  layer  becomes  either  laminar  or 
transitional  depending  on  rotation  number  and  buoyancy  forces.  In  radial 
outward  flow  for  the  cases  examined,  the  heat  transfer  becomes  laminar  for 
10^®  <  Rajj  <  10^^  and  can  be  predicted  with  a  flat  plate  laminar  heat 
transfer  correlation.  Outside  of  this  range,  the  boundary  layer  becomes 
transitional  and  is  most  difficult  to  accurately  correlate. 


8.3  Comparison  of  Results  with  Other  Experiments 


Many  authors  have  examined  the  effects  of  rotation  from  an  Isothermal, 
hydrodynamic  point  of  view.  However,  few  have  experimentally  analyzed  the 
effects  of  rotation  on  internal  heat  transfer  and  fewer  still  have  measured 
the  effects  of  centrifugal  buoyancy  and  Coriolis  forces  on  turbulent  heat 
transfer  for  conditions  similar  to  those  found  in  today's  advanced  gas  turbine 
blades.  The  following  sections  compare  the  heat  transfer  data  of  this  program 
to  the  results  of  several  published  experiments. 
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Figure  57.-  Leading  Side  Heat  Transfer  with  Rotation  Compared  to  Kays  (1966) 
Laminar  Heat  Transfer:  (Re  =  12,500,  R/d  =  49,  Q!=  0). 
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Figure  58.-  Leading  Side  Heat  Transfer  with  Rotation  Compared  to  Kays  (1966) 
Laminar  Heat  Transfer:  (Re  =  50,000,  R/d  =  49,  a=  0). 
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Figure  59.-  Present  Data  and  Johnston's  Flow  Map  (Johnston,  1970)  Showing 
Laminarization  on  the  Low  Pressure  Surface. 


8.3.1  Comparison  with  Results  of  Iskakov  (1985) 

The  upturn  in  heat  transfer  rates  measured  on  the  low  pressure  surface  with 
increasing  rotation  is  a  phenomena  not  reported  until  recently.  Only  Iskakov 
(1985)  is  known  to  have  documented  this  trend.  In  a  rotating  circular  duct 
experiment,  he  acquired  local  heat  transfer  coefficients  on  both  the  leading 
and  trailing  surfaces  for  radially  outward  flow  for  several  Reynolds  numbers. 
At  Reynolds  number  Re=32000,  Iskakov  reported  both  a  decrease  in  the  low 
pressure  surface  heat  transfer  to  about  85%  of  the  stationary  level  and  a 
significant  upturn  as  rotation  increased.  He  showed  this  trend  is  strongly 
Reynolds  number  dependent.  As  expected,  on  the  other  side  of  the  circular 
duct,  the  high  pressure  side,  Iskakov  measured  increasing  heat  transfer  rates. 
These  results  are  presented  in  Figure  60.  Iskakov 's  result  are  compared  to 
data  of  this  program  in  Figure  61.  His  correlations  are  plotted  with  data  for 
Re=25000  at  six  rotation  rates.  The  heat  transfer  coefficients  have  been 
averaged  over  the  passage  length  for  the  high  and  low  pressure  sides  of  the 
model  for  radially  outward  flow.  Recall  that  Iskakov 's  duct  is  circular  and 
the  present  duct  is  square. 
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Figure  61.-  Comparison  of  Iskakov's  Correlation  with  Averaged  Heat  Transfer 
Data  of  this  Experiment:  Re  =  25,000  (Iskakov  and  Trushin,  1985) 


On  the  high  pressure  side  of  both  models,  the  enhancement  in  heat  transfer 
with  increasing  rotation  agrees  up  to  a  rotation  number  of  Ro=0.30.  On  the  low 
pressure  surface  both  exhibit  a  drop  and  then  an  Increase  in  heat  transfer, 
but  the  levels  remain  very  different. 

This  suggests  that  a  square  duct  may  be  more  susceptible  than  a  circular  tube 
to  decreases  in  heat  transfer  on  the  low  pressure  side.  This  is  believed  to 
result  from  the  Coriolis-induced  stabilization  of  the  wall  shear  layer. 

Because  there  are  corners  in  the  square  duct,  the  secondary  flows  are  less 
likely  to  wash-off  this  stabilized  fluid  layer.  It  is  suspected  that  the 
secondary  circulations  migrate  toward  the  high  pressure  surface  leaving  a 
thickened  and  dampened  boundary  layer  on  the  low  pressure  surface.  This  would 
explain  why  the  square  duct  experiences  such  a  large  drop  in  heat  transfer  on 
this  surface;  the  level  approaching  that  of  laminar  flat  plate  heat  transfer. 

It  is  interesting  to  point  out  that  both  experiments  report  an  upturn  in  low 
pressure  surface  heat  transfer  at  a  rotation  number  near  Ro  =  0.20.  This 
upturn  in  heat  transfer  on  the  low  pressure  side  of  the  rotating  duct  may  be  a 
flow  separation  phenomena  (see  Iskakov  and  Trushin,  1983  and  1985). 

8.3.2  Comparison  with  Results  of  Morris  (1981  ) 

Comparisons  with  the  data  and  correlation  of  Morris  illustrate  one  of  the  most 
important  realizations  of  the  program:  Extrapolation  of  correlations  outside 
the  test  range  of  the  significant  nondimens ional  parameters  is  very  risky. 

Figure  62  compares  Morris'  correlation  to  data  from  the  first  straight  passage 
in  the  present  model.  Morris'  experiment  consisted  of  a  rotating  circular  tube 
with  radially  outward  flow.  This  was  a  constant  heat  flux  experiment  where  a 
circular  duct  was  instrumented  with  thermocouples  along  the  leading  surface. 
These  results  are  then  compared  to  those  from  the  first  passage  at  X/d  =  12.4 
(radially  outward  flow)  on  the  leading  surface  of  the  NASA  model. 

The  solid  lines  of  Figure  62  indicate  Morris'  correlation  based  on  data,  while 
the  dashed  line  represent  this  correlation  extrapolated.  The  heat  transfer 
results  obtained  on  the  leading  surface  at  the  low  rotation  number  of  0.06  are 
within  20%  of  Morris'  correlation.  This  agreement  occurs  within  Morris'  experi 
mental  range  which  is  indicated  in  the  figure.  For  values  of  the  buoyancy  para 
meter  and  rotation  number  that  fall  outside  the  Morris  data  range,  the  present 
data  show  increases  with  rotation  niimber  which  are  in  general  agreement  with 
the  Morris  correlation.  However,  the  present  data  show  increases  in  heat  trans 
fer  with  increasing  density  ratio  or  buoyancy  centripetal  parameter,  whereas, 
the  Morris  correlation  would  predict  a  decrease  in  heat  transfer  with 
increasing  density  ratio. 

8.3.3  Comparison  with  Results  of  Eckert  (1954') 

This  program  has  shown  rotational  buoyancy  forces  strongly  influence  turbulent 
heat  transfer  in  rotating  smooth  passages  for  conditions  found  in  gas  turbine 
blades.  Variations  in  heat  transfer  caused  by  rotation  have  been  shown  to  be 
less  for  radial  inward  flowing  passages  then  for  radial  outward  flowing 
passages.  This  influence  of  radial  flow  direction  is  clarified  when  the 
results  of  Eckert  and  Diaguilla's  (1954)  work  on  buoyancy  effects  in  forced, 
free  and  mixed  convection  are  examined. 
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Figure  62.-  Comparisons  of  Leading  Surface  Rotating  Heat  Transfer  Results  for 
X/d  =  12.4  with  Morris  (1981)  Correlation: 

Nu  =  0.022  (Ra/Re2)-0*186Ro0.33Re0-8  or 

Nu/NUoo=  1.35  [(AP/P)*(R/d)*(n/d/V)2]-0.186(nd/v)0.33 


Eckert  conducted  three  different  experiments  in  a  large  (stationary)  vertical 
circular  tube  where  he  examined  mixed,  free  and  forced-convection  heat 
transfer.  A  diagram  of  the  test  apparatus  is  shown  with  the  results  of  his 
experiments.  These  are  plotted  as  Nu^  versus  Gr^  Pr  in  Figures  63  and  64. 

Each  of  these  plots  contains  a  free-convection  limit  which  is  a  mean  data  fit 
through  the  results  of  Eckert's  free-convection  experiment. 

Figure  63  compares  the  present  data  (symbols)  for  the  outward  flowing  passage 
to  Eckert's  curve-fitted  (solid)  lines  for  his  counterflow  tests.  In  this  case 
Eckert  looked  at  the  condition  where  the  wall  thermal  boundary  layer  under 
gravity-induced  buoyancy  forces  moves  in  a  direction  opposite  to  the  mainstream 
flow.  This  is  analogous  to  the  rotating  case  of  the  outward  flowing  passage 
with  heated  walls  where  the  buoyancy  force  is  rotation- induced.  Here  the  more 
dense  colder  fluid,  away  from  the  wall,  accelerates  outward  from  the  axis  of 
rotation,  while  the  hotter  wall  fluid  tends  towards  the  axis  of  rotation.  All 
heat  transfer  coefficients  presented  are  averaged  circumferentially  around  the 
perimeter  of  the  respective  passage. 
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Figure  63.- 


Circxjinferentially  Averaged  Heat  Transfer  Coefficients  for  Outflow 
Passage  Compared  with  Eckert's  Counter flow  Results  (Eckert  and 
Diaguilla,  1954). 


STEAM 

JACKET 


ECKERT'S  MIXED,  FREE,  AND 
FORCED-CONVECTION  TUBE 


Figure  64.-  Circumferentially  Averaged  Heat  Transfer  Coefficients  for  Inflow 
Passage  Compared  with  Eckert's  Parallel  Flow  Results  (Eckert  and 
Diaguilla,  1954). 
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The  present  data  plotted  for  various  Reynolds  ntimbers  based  on  the  distance  X, 
measured  from  the  beginning  of  the  thermal  boundary  layer,  shows  very 
encouraging  trends.  Moving  to  the  left  most  side  of  Figure  63,  the  lines  of 
constant  Re^  should  eventually  become  flat  as  the  heat  transfer  becomes 
independent  of  Grasshof  number  and  can  be  correlated  by  an  equation  of  the  form 
Nujj  =  CRe^Pr^.  Eckert  noted  for  this  case  that  the  data  did  not  yet 
become  flat,  and  therefore,  the  heat  transfer  at  these  lower  Grasshof  niunber 
remains  dependent  on  buoyancy.  Moving  to  the  right  side  of  the  figure,  for 
large  Grasshof  numbers,  heat  transfer  rates  approach  the  free  convection  limit 
established  by  Eckert. 

The  present  data,  when  averaged  circumferentially  around  the  duct  perimeter, 
agrees  remarkably  well  with  Eckert's  results.  The  increasing  slope  through  the 
data  indicates  the  flow  is  in  the  mixed  flow  regime  and  buoyancy  influences 
the  heat  transfer  process.  For  large  values  of  Gr^,  heat  transfer  rates 
approach  the  free  convection  limit  as  established  by  Eckert.  The  data  also 
indicates  a  higher  free  convection  limit  at  the  lower  values  of  Gr^.  Eckert 
also  noticed  this  trend  when  comparing  the  extreme  forced-convection  data  at 
high  Grashof  nrombers  with  the  free  convection  data.  For  low  values  of  Gr^, 
forced  convection  dominates  the  heat  transfer.  As  can  be  seen  in  the  figure, 
the  forced  convection  limit  was  not  reached  with  the  present  set  of 
experiments . 

Figure  64  compares  the  present  results  for  the  inward  flowing  passage  with 
Eckert's  parallel  flow  results.  Here,  Eckert  reversed  the  mainstream  flow 
direction  producing  a  thermal  boundary  layer  that  moved  in  the  same  direction 
as  the  mainstream.  Again,  the  two  experiments  compare  favorably.  The  present 
data  falls  within  the  mixed  regime  defined  by  Eckert. 

This  mixed  region  exists  between  the  forced-convection  and  free-convection 
limits  indicated  on  Figure  64.  It  is  important  to  notice  that  for  this 
parallel  flow  case,  Eckert's  data  established  the  pure  forced-convection 
regime  where  Nu^  is  independent  of  Gr^.  As  previously  noted,  the 
counterflow  experiments  did  not  show  this  flat  region  (see  Figure  63).  This 
suggests  that  in  the  parallel  flow  case,  buoyancy  influences  heat  transfer 
over  a  smaller  range.  Thus,  for  radially  outward  flowing  passages,  heat 
transfer  will  more  likely  be  affected  over  a  wider  operating  range  by 
rotation- induced  buoyancy  than  for  radially  inward  flowing  passages. 

It  has  been  shown  that  thermal  boundary  layer  buoyancy  acting  in  the  S2une 
direction  as  the  turbulent  mainstream  inhibits  heat  transfer  in  the  mixed 
convection  region.  Both  Eckert's  experiments  in  the  stationary  mode  and  the 
rotating  results  presented  here  exhibit  a  local  decrease  in  Nu^  when  Gr^Pr 
falls  in  the  mixed  range  shown  in  Figure  64.  This  local  decrease  is  especially 
significant  when  compared  with  counterflow  case.  Vflien  Nusselt  numbers  for 
similar  flow  conditions  for  the  radially  outward  and  inward  flowing  case  are 
compared,  the  counterflow  (radially  outflow)  heat  transfer  is  almost  70% 
greater  than  the  corresponding  parallel  flow  heat  transfer  (e.g.  Nu^  =  875 
compared  to  520).  Eckert  noted  that  the  counterflow  situation  could  result  in 
heat  transfer  levels  as  much  as  twice  those  in  the  parallel  flow  case.  This  is 
especially  significant  when  these  results  are  circumferential  averages  of  heat 
transfer  around  the  perimeter  of  the  coolant  passage.  As  shown  in  the  previous 
discussion,  local  heat  transfer  can  be  significantly  higher  and  lower  than  the 
circumferential  averages. 
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Because  of  the  absence  of  detailed  boundary  layer  measurements  for  these 
experiments,  it  is  not  clear  what  turbulence  structure  and/or  mean  flow 
changes  are  responsible  for  the  observed  heat  transfer  behavior.  Clearly, 
measurements  of  this  nature  will  be  required  to  further  the  physical 
understanding  of  this  class  of  shear  flows  and  thereby  permit  the  development 
of  more  accurate  predictive  schemes. 

Several  conclusion  regarding  buoyancy  can  be  reached: 

1.  Rotating  coolant  passage  heat  transfer,  averaged  circumferentially  for  a 
square  smooth  channel,  compares  favorably  with  Eckert's  work  on 
simultaneously  occurring  free  and  forced-convection  in  a  stationary 
vertical  circular  tube. 

A.  The  radial  outward  flowing  passage  with  rotation  correlates  with  the 
stationary  counterflow  case. 

B.  The  radial  inward  flowing  passage  with  rotation  correlates  with  the 
stationary  parallel  flow  case. 

2.  The  large  average  heat  transfer  variations  seen  between  the  outward  and 
inward  flowing  passages  are  attributed  to  buoyancy  forces. 

3.  The  mixed  convection  regime,  where  free  and  forced-convection  occur 
simultaneously,  is  larger  for  radially  outward  flowing  passages  (opposed 
flow)  than  for  inward  flow  (parallel  flow).  As  a  consequence,  heat 
transfer  is  more  affected  by  rotation  in  the  outward  flowing  passage. 
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9.0  CONCLUSIONS 


This  progam  has  resulted  in  an  extensive  body  of  experimental  data  from  heat 
transfer  experiments  in  a  rotating  square  passage  with  smooth  walls.  It  is 
believed  that  the  large  range  of  test  parameters  makes  this  data  set  unique. 
The  extensive  data  base  aided  greatly  in  the  data  analysis  and  correlation  and 
in  developing  physical  models  for  the  complex  heat  transfer  characteristics.  A 
summary  of  the  major  program  conclusions  concerning  the  separate  effects  of 
forced  convection,  Coriolis,  buoyancy  and  flow  direction  on  heat  transfer  is 
presented  in  this  section.  A  more  detailed  discussion  of  these  conclusions  and 
the  authors'  speculations  regarding  physical  models  for  the  cause  and  effect 
relationships,  can  be  found  in  Sections  5.0,  6.0,  7.0  and  8.0  of  this  report. 

For  the  first  passage  with  radially  outward  flow: 

1.  The  streamwise  distribution  of  heat  transfer  in  the  first  passage  from 
the  stationary  baseline  experiment  is  similar  to  that  for  developing 
flow  in  the  entrance  region  of  a  passage  and  is  in  good  agreement  with 
previous  investigators'  results. 

2.  Heat  transfer  is  strongly  affected  by  rotation,  causing  increases  in 
heat  transfer  up  to  3.5  times  fully  developed,  smooth  tube  levels  on  the 
trailing  surfaces  and  decreases  to  40%  of  fully  developed,  smooth  tube 
levels  on  the  leading  surfaces. 

3.  The  decreases  in  heat  transfer  on  the  leading  surfaces  with  increases  in 
rotation  number  are  attributed  to  the  combined  effects  of  stabilization 
of  the  near-wall  flow  and  cross-stream  flows  which  cause  heated,  near 
wall  fluid  from  the  trailing  and  sidewall  surfaces  to  accumulate  near 
the  leading  side  of  the  coolant  passage. 

4.  The  increases  in  heat  transfer  in  the  latter  half  of  the  coolant  passage 
on  the  leading  side  and  the  Increases  on  the  trailing  side  are 
attributed  to  the  effects  of  the  large  scale  development  of  Coriolis 
generated  secondary  flows. 

5.  Local  heat  transfer  increases  (by  as  much  as  200%)  with  increases  in 
density  ratio  over  most  of  the  passage  surface  area. 

6.  Heat  transfer  decreases  (by  as  much  as  20%)  with  increases  in  density 
ratio  on  the  leading  side  of  the  passage  near  the  inlet.  These  decreases 
are  believed  to  be  governed  by  both  the  interaction  of  the  near-wall 
flow  stability  and  the  buoyancy  effects. 

7.  Heat  transfer  results  obtained  with  different  combinations  of  model 
radius  ratio  and  coolant  density  were  well  correlated  by  the  buoyancy 
parameter  (AP/P ) (R/d)  (nd/V)2  for  flow  regimes  dominated  by  buoyant 
forces . 
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8. 


Increases  in  heat  transfer  ratio  on  the  trailing  surfaces  were  60% 
greater  than  increases  on  the  leading  surfaces  for  the  same  increase  in 
density  ratio  for  a  given  rotation  number.  This  difference  in  heat 
transfer  increase  suggests  that  the  interaction  of  the  Coriolis  and 
buoyancy  effects  is  different  for  the  flow  near  leading  and  trailing 
surfaces. 

9.  The  buoyancy  parameter  correlates  the  heat  transfer  ratio  results  from 
the  trailing  side  and  for  the  most  downstream  location  of  the  leading 
side.  The  results  were  not  correlated  by  the  buoyancy  parameter  near  the 
inlet  on  the  leading  surface.  The  lack  of  correlation  was  attributed  to 
a  complex  interaction  of  stabilization  of  the  near  wall  flow,  buoyancy 
and  Coriolis  effects. 

The  following  observations  and  conclusions  were  reached  concerning  flow 
direction: 

1.  The  effects  of  varying  Reynolds  nximber  on  heat  transfer  in  the  first  two 
passages  are  reasonably  well  correlated  by  the  heat  transfer  ratio 
Nu/Nu  ,  where  Nu  is  the  correlation  for  fully  developed,  turbulent  flow 
in  a  stationary  rectangular  passage. 

2.  Density  ratio  and  rotation  nvimber  were  found  to  cause  large  changes  in 
heat  transfer  ratio  for  radially  outward  flow  and  relatively  small 
changes  for  radially  inward  flow. 

3.  The  heat  transfer  ratio  for  the  low  pressure  surfaces  of  the  coolant 
passages  was  found  to  be  primarily  a  function  of  the  buoyancy  parameter 
regardless  of  flow  direction. 

4.  The  heat  transfer  ratio  on  the  high  pressure  surfaces  was  significantly 
affected  by  flow  direction.  The  heat  transfer  was  found  to  be  a  strong 
function  of  a  buoyancy  parameter  for  the  high  pressure  surfaces  for 
radially  outward  flow.  Whereas,  the  heat  transfer  was  relatively 
unaffected  by  a  buoyancy  parameter  for  the  radially  inward  flowing  high 
pressure  surface. 

5.  Increasing  the  density  ratio  generally  caused  an  Increase  in  heat 
transfer.  However,  the  Increase  in  heat  transfer  for  the  inward  flowing 
passage  was  considerably  less  than  that  for  outward  flow. 

6.  Circiimferentially  averaged  heat  transfer  results  compared  favorably  with 
a  previous  stationary  parallel  and  counterflow  mixed-convection  heat 
transfer  experiment  by  Eckert  et  al.  (1953).  Heat  transfer  from  the 
radially  outward  flowing  passages  compared  with  Eckert’s  counterflow 
case  while  the  radially  inward  flowing  results  compared  with  Eckert's 
parallel  flow  case. 

7.  The  effect  of  model  orientation  angle  ( a  =  45®)  caused  heat  transfer  to 
be  symmetric  about  a  diagonal  in  the  direction  of  rotation  across  the 
passage. 
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10.0  APPENDIX 


10.1  ANALYSIS  OF  EQUATIONS  OF  MOTION  (Suo,  1980) 

Consider  a  fluid  flowing  in  a  pipe  which  is  rotating  as  shown  in  Figure  65. 


Figure  65.-  Rotating  Pipe  Coordinates. 


The  equation  of  motion  in  a  rotating  frame  of  reference  is 


P  -  Vp  +  /fV(V.V)  +  V-  p  [(2^  X  P)  +  (£I  X  S?  X  r)] 

^  Coriolis  Centrifugal 


(10.1) 


The  last  two  terms  represent  the  Coriolis  and  centrifugal  acceleration  terms, 
respectively.  If  Eq.  (10.1)  is  made  dimensionless  in  the  normal  manner. 
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Then  Equation  (10.1)  becomes 
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where  Re  is  (V^dP^/M)  and  (dfi/V^)  is  the  rotation  number,  (i.e.  inverse 
of  Rossby  number) .  In  order  to  determine  the  important  variables  in  the 
problem  of  concern,  Eq.  (10.3)  is  further  modified  by  the  definition  of 
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Equation  (10.3)  becomes 
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The  relative  importance  of  each  of  the  terms  on  the  right  can  then  be 
examined.  First  consider  the  case  where  (AP/P^)  is  very  large.  Then  all 
terms  are  important  and  the  solution  becomes  very  difficult.  An  example  of 
such  a  flow  is  the  flow  in  a  rotating  pipe  where  boiling  occurs.  In  such  a 
flow  the  liquid  has  been  observed  to  stratify  on  the  side  of  the  pipe  away 
from  the  direction  of  rotation.  This  flow,  while  interesting  in  its  own  right, 
is  not  of  interest  for  the  problem  of  air-cooled  turbine  blades. 


A  second  extreme  is  whenAP/Pm  =  0.  Then  Eq.  (10.6)  reduces  to 


D7' 

Dt' 


If  we  define  a  reduced  pressure  gradient  as, 
Vp*  =  Vp-  P„rQ^t  =  v[p-1/2  P„  (^f] 


Then  Equation  (10.7)  becomes 


(10.7) 


(10.8) 


(10.9) 


The  reduced  pressure  gradient  Vp*  is  the  pressure  gradient  minus  the  pressure 
gradient  caused  by  centrifugal  acceleration.  Thus,  it  represents  the  pressure 
gradient  driving  the  flow,  because,  if  only  the  centrifugal  gradient  was 
present,  there  would  be  no  motion  except  solid  body  rotation.  The  only  term  to 
make  this  equation  different  from  nonrotating  flows  is  the  last  one,  which 
represents  the  effect  of  Coriolis  acceleration. 


The  flows  of  most  interest  for  the  current  problem  are  not  at  either  of  the 
extremes,  but  lie  in  between.  Consider  the  case  when  P/Pjn  <<  1  hut  is  not 
equal  to  zero.  This  allow  the  Bousslnesq  approximation  to  be  made  where  the 
variation  in  density  is  only  considered  in  the  driving  terms,  here  the 
rotational  terms 
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Again  using  the  definition  of  reduced  pressure  gradient,  Eq  (10.10)  becomes 
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The  first  two  terms  within  the  brackets  are  the  Coriolis  acceleration  terms. 
Clearly  the  second  term  is  always  very  much  smaller  than  the  first  term  and 
can  be  neglected  with  respect  to  the  first  term.  The  third  term  is  the 
centrifugal  acceleration  term;  it  can  approach  the  value  of  the  first  term 
(with  rD/Vm>>l)  and  thus  cannot  be  neglected  in  general.  Equation  (10.11) 
then  becomes 
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The  term  on  the  far  right  represents  the  effect  of  centrifugal  buoyancy.  When 
(r^Vv^)  approaches  zero,  the  effect  of  centrifugal  buoyancy  should 
become  negligible  and  Equation  (10.12)  reduces  to  Equation  (10.9)  above. 
Furthermore,  it  is  clear  that  in  addition  to  Reynolds  number,  the  terms 
(dfl/Vjn)  and  AP/Pm(rn/Vni)  should  specify  the  flow  field. 

10.2  Model  Inlet  Flow  Characterization 

Velocity  measurements  were  obtained  at  the  exit  of  the  inlet  nozzle  system 
(inlet  of  the  first  passage)  and  at  the  end  of  a  square  chaimel  mounted  on  the 
inlet  nozzle  (Figure  8).  The  velocity  measurements  were  obtained  by  traversing 
a  hot  wire  probe  across  the  passage  from  the  leading  to  trailing  side  and  from 
the  inside  of  the  passage  to  the  outside.  Average  and  RMS  voltage  output  of 
the  linearized  hot  wire  signal  were  used  to  calculate  the  local  mean  velocity 
and  local  turbulence  intensity. 

Probe  traverse  locations  and  nomenclature  are  shown  in  Figure  66.  Mean 
velocity  and  local  turbulence  profiles  from  the  measurements  are  shown  in 
Figures  67  through  70.  The  mean  velocity  results  for  the  flow  at  the  nozzle 
exit  (entrance  to  the  first  passage)  (Figures  67a  and  67b  top)  show  the 
velocity  profiles  have  an  approximately  parabolic  shape.  The  mean  velocity 
surface  and  contour  plots  for  the  nozzle  exit  (Figure  68a)  show  the  core  flow 
to  be  approximately  axisymmetric  about  the  passage  centerline  and  slightly 
skewed  toward  the  outside  and  toward  the  trailing  side  of  the  passage. 
Turbulence  profiles  at  the  nozzle  exit  (Figures  67a  and  67b  bottom  and  68b) 
show  increased  local  turbulence  intensity  near  the  edges  of  the  nozzle  flow 
with  a  centerline  turbulence  of  about  4%. 

Mean  velocity  profiles  at  the  end  of  the  channel  (Figures  69a  and  69b  top)  are 
relatively  flat  and  typical  of  fully  developed  turbulent  flow.  Surface  and 
contour  results  (Figure  70a)  indicate  that  the  mean  flow  is  slightly  skewed 
toward  the  leading  side  of  the  passage.  Turbulence  profiles  at  the  channel 
exit  (Figures  69a  and  69b  bottom,  and  70b)  show  a  decreased  level  of 
centerline  turbulence  intensity,  down  from  4%  at  the  nozzle  to  about  3%  at  the 
channel  exit. 
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TRAILING  SIDE 


LOCAL  TURBULENCE  INTENSITY,  u/Uj  MEAN  VELOCITY  RATIO,  U/U 


a)  MEAN  VELOCITY  AND  LOCAL 
TURBULENCE  IN  Z  DIRECTION 
□y  =  0.05  Ay  =  0.15  oy  =  0.25 
xy  =  0.35  +y  =  0.45 


b)  MEAN  VELOCITY  AND  LOCAL 
TURBULENCE  IN  Y  DIRECTION 
Dz  =  0.05  A2  =  0.15  Oz  =  0.25 
xz  =  0.35  +z  =  0.45 
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Figure  67.-  Nozzle  Exit  Flow  Profiles  in  Z  and  Y  Directions 
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<  LEADING  SIDE, 


a)  MEAN  VELOCITY  SURFACE 
AND  CONTOUR  ,  Uj/U 


b)  LOCAL  TURBULENCE  SURFACE 
AND  CONTOUR  ,  U|7U| 


OUTWARD  SIDE 


Figure  68.-  Variation  in  Velocity  Surfaces  and  Contours  for  Nozzle  Exit 
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LOCAL  TURBULENCE  INTENSITY,  Ui/Uj  MEAN  VELOCITY  RATIO ,  U/U 


a)  MEAN  VELOCITY  AND  LOCAL  b)  MEAN  VELOCITY  AND  LOCAL 
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Figure  69.-  Channel  Exit  Flow  Profiles  in  Z  and  Y  Directions 


a)  MEAN  VELOCITY  SUR^FACE  b)  LOCAL  TURBULENCE  SURFACE 

AND  CONTOUR  ,  U/U  AND  CONTOUR  ,  u/Uj 


Figure  70.-  Variation  in  Velocity  Surfaces  and  Contours  for  Channel  Exit  Flow. 


10.3  Error  Analysis,  Repeatability  and  Symmetry  of  the  Model  Heat  Transfer 
Results 

An  analysis  to  determine  the  estimated  error  in  measured  heat  transfer  was 
conducted  using  the  method  described  in  Abernathy  and  Thompson  (1973).  The 
error  analysis  was  a  summation  of  the  estimated  inaccuracies  in  the  data  used 
to  calculate  the  heat  transfer  coefficient  and  Nusselt  number  as  shown  below. 
When  comparative  testing  is  done  (as  in  this  program),  Abernathy  and  Thompson 
suggest  that  biasing  errors  may  be  ignored.  Therefore,  only  the  precision 
errors  were  determined. 

n 

(ANu)"  =  Y  lldNu/df,)"  (Af,  )2]  (10.13) 

/  =  / 

where  f^  are  the  parameters  used  to  determine  the  Nusselt  number 

To  do  the  error  analysis,  assumptions  of  the  estimated  inaccuracies  had  to  be 
made.  Generally,  where  information  was  read  by  eye  from  a  device  (i.e.,  flow¬ 
meter,  pressure  guage,  etc.),  the  estimated  accuracy  was  assumed  to  be  1/4  of 
the  smallest  increment.  Voltages  read  by  the  acquisition  system  were  assumed 
to  be  accurate  to  the  least  significant  digit,  except  for  the  thermocouple 
emfs  which  were  assumed  to  be  accurate  to  only  25  micro-volts.  A  summary  of 
the  estimated  inaccuracies  follows: 


■ameter 

Estimated  Error  <’+') 

Method 

VD 

0.000025  Volt 

T.C.  wire  repeatability 

VC 

0.000025  Volt 

T.C.  wire  repeatability 

VB 

0.0001  Volt 

least  significant  measurement 

E 

0.001  Volt 

least  significant  measurement 

I 

0.001  Volt 

least  significant  measurement 

QALOSS 

(0.05)  QALOSS 

5%  accurate  estimate 

FMA 

0.25 

1/4  Increase  on  flow  meter 

CFMA 

(0.005)CFMA 

0.5%  flow  meter  reading 

BFMA 

(0.005)BFMA 

0.5%  flow  meter  reading 

PMin 

0.25 

1/4  increase 

TMin 

0.5 

1/2  resolution 

h 

(0.02)h  at  the  inlet 

2% 

(0.06)h  at  the  exit 

6% 

Note:  h  was  only  used  to  estimate  the  error  of  determining  the  effective 
heat  transfer  surface  area. 

This  analysis,  applied  to  the  case  where  the  stationary  test  point  Reynolds 
Number  is  25,000,  showed  the  error  of  determing  heat  transfer  at  the  entrance 
of  the  model  was  +  2%  and  the  error  at  the  exit  was  ±  7%  of  the  local  heat 
transfer  rate.  Of  the  estimated  error,  over  90%  was  due  to  the  error  in 
determining  the  temperature.  Nine  percent  of  the  estimated  error  was  due  to 
the  error  in  determining  the  backloss. 

A  series  of  tests  were  conducted  to  determine  the  accuracy,  symmetry  and 
repeatability  of  the  rotating  heat  transfer  experiments.  The  data  acquisition 
procedure.  Instrumentation,  model  calibration  and  heat  loss  aspects  of  the 
accuracy  problem  were  addressed  in  Section  3.  The  results  from  the  tests  to 
measure  repeatability  and  model  symmetry  are  reported  here. 
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10.3.1  Repeatability 

The  repeatability  of  the  heat  transfer  experiments  for  the  baseline  rotating 
flow  condition  was  determined  by  comparing  the  heat  transfer  coefficients 
obtained  from  UTRC  Run  22.12  (Test  No.  4)  and  UTRC  Run  25.08.  The  "repeat" 
test  was  conducted  after  Test  Nos.  9  and  10  which  had  lower  and  higher  wall 
temperatures  than  the  baseline  flow  condition.  The  ratio  of  the  heat  transfer 
coefficients  from  the  baseline  flow  condtion  (Test  4)  to  that  from  the  repeat 
test  is  shown  in  Table  IV  for  each  of  the  heater  elements.  The  location  for 
each  heater  element  is  identified  in  Figure  71.  A  comparison  of  the  elements 
from  the  first  leg  of  the  model,  i.e.  elements  1-4,  19-22,  33-36  and  49-52, 
shows  that  the  variation  of  the  "repeat"  results  from  the  baseline  results  are 
1  or  2%,  except  for  one  element.  The  variation  for  the  repeat  results  in  the 
second  leg  of  the  model  is  increased  to  4%  and  the  variation  in  the  third  leg 
is  up  to  9%.  The  increase  in  the  variation  of  the  results  as  the  distance  from 
the  inlet  Increased  was  attributed  to  the  uncertainty  in  the  bulk  temperature 
which  was  calculated  as  part  of  the  data  reduction  procedure.  The  repeatability 
was  considered  very  good,  considering  the  large  decrease  in  heat  transfer  due 
to  rotation  on  the  leading  surfaces  of  the  first  and  third  passages. 

10.3.2  Low  Rotation 


Tests  were  also  conducted  at  a  low  rotational  rate  (±15  rpn./ .  The  results  from 
the  tests  at  15  rpm  (Tests  101  and  108)  are  compared  with  results  from  the 
stationary  test  (lA)  at  the  baseline  operating  flow  condition  as  shown  in 
Figure  72.  Note  that  for  this  flow  condition  with  dimensionless  rotation 
number  (Ro)  equals  -0.006,  0,  and  +0.006,  small  effects  of  rotation  are  noted. 
However,  the  results  with  rotation  are  either  very  close  (i.e.  5%)  or  bracket 
the  zero  rotation  results.  The  conclusion  from  this  comparison  is  that 
acquisition  of  data  through  the  slipring  does  not  significantly  affect  the 
heat  transfer  results. 

10.3.3  Effect  of  Rotation  Direction 

Data  was  also  acquired  with  the  model  rotating  in  both  the  positive  and 
negative  directions  at  the  baseline  operating  conditions.  These  tests  were  con¬ 
ducted  to  determine  if  small  asymmetries  in  the  construction  of  the  elements  or 
the  assembly  of  the  model  were  affecting  the  heat  transfer  results.  The  results 
from  Test  Nos.  4  (+550  rpm)  and  109  (-550  rpm)  are  presented  in  Figure  73.  The 
general  results  from  both  tests  are  identical.  However,  there  is  as  much  as  a 
±10%  variation  in  the  heat  transfer  coefficient  obtained  from  different  ele¬ 
ments  in  the  s^^ne  relative  location.  For  example,  element  35  (upper  left  por¬ 
tion  of  Figure  73)  at  the  C  streamwise  location  on  the  leading  surface  has 
ratios  of  Nu/Nu^q  =  0.45  and  0.53,  respectively.  The  conclusions  from  this 
comparison  of  results  from  the  ±  550  rpm  test  conditions  are  (1)  that  there 
are  no  major  asymmetries  in  the  model,  and  (2)  that  the  principal  results  from 
the  test  are  not  sensitive  to  model  asymmetries. 
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TABLE  IV 

SYMMETRY  OF  MODEL  AND  REPEATABILITY 

OF  HEAT  TRANSFER 

RESULTS 

ROTATING  HEAT 

TRANSFER  EXPERIMENTS 

WITH  SMOOTH  WALL 

MODEL 

(Run  No:  22:12  and 

25.08) 

Dimensionless  Data 

Element 

No.  ~X/d 

Ro 

Nn  r22.12')/Nu  (25 

1 

1 

0.2410 

0.99 

2 

5 

0.2386 

0.99 

3 

9 

0.2365 

0.98 

4 

12 

0.2345 

0.98 

5 

14 

0.3106 

0.97 

6 

15 

0.3106 

0.99 

7 

16 

0.3079 

0.98 

8 

17 

0.3079 

0.99 

9 

2 

0.2286 

0.99 

10 

6 

0.2274 

0.99 

11 

10 

0.2266 

0.97 

12 

2 

0.2245 

0.98 

13 

6 

0.2237 

0.97 

14 

10 

0.2230 

0.98 

15 

11 

0.2973 

0.95 

16 

12 

0.2973 

0.94 

17 

14 

0.2964 

1.03 

18 

15 

0.2964 

0.95 

19 

1 

0.2410 

1.00 

20 

5 

0.2386 

0.98 

21 

9 

0.2365 

0.98 

22 

12 

0.2345 

0.99 

23 

2 

0.2286 

0.99 

24 

6 

0.2274 

0.98 

25 

10 

0.2266 

0.98 

26 

11 

0.3021 

0.99 

27 

12 

0.3021 

1.00 

28 

14 

0.3012 

0.99 

29 

15 

0.3012 

0.97 

30 

2 

0.2245 

0.98 

31 

6 

0.2237 

0.99 

32 

10 

0.2230 

0.99 

33 

1 

0.2410 

1.00 

34 

5 

0.2386 

1.01 

35 

9 

0.2365 

1.01 

36 

12 

0.2345 

0,95 

37 

14 

0.3106 

0.98 

38 

17 

0.3079 

0.98 

39 

2 

0.2286 

0.96 

40 

6 

0.2274 

0.97 

41 

10 

0.2266 

0.99 
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TABLE  IV  (Cont’d) 


Dimensionless  Data 


Element  No. 

~X/d 

Ro 

Nu  (22.12)/m  C25.081 

42 

12 

0.3021 

0.97 

43 

14 

0.3012 

0.98 

44 

2 

0.2245 

0.96 

45 

6 

0.2237 

0.94 

46 

10 

0.2230 

0.92 

47 

12 

0.2973 

0.91 

48 

14 

0.2964 

0.96 

49 

1 

0.2410 

1.00 

50 

5 

0.2386 

0.98 

51 

9 

0.2365 

0.98 

52 

12 

0.2345 

0.98 

53 

14 

0.3106 

0.99 

54 

17 

0.3079 

0.98 

55 

2 

0.2286 

0.97 

56 

6 

0.2274 

0.96 

57 

10 

0.2266 

0.96 

58 

12 

0.3021 

0.95 

59 

14 

0.3012 

0.97 

60 

2 

0.2245 

0.97 

61 

6 

0.2237 

0.98 

62 

10 

0.2230 

0.97 

63 

12 

0.2973 

0.94 

64 

14 

0.2964 

0.93 
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HEAT  TRANSFER  RATIO,  Nu/Nu 


Re  ^25,000 
Ro^.24 
RPM  =  550 


AT  =  44.4  C®(80°F) 
R/d  =  49 


SYMBOL 

TEST  NO 

UTRC  RUN  NO 

o 

▲ 

4 

Repeat 

22.12 

25.8 

LEADING  (  +  Q) 
SEGMENTS  33-48 


TRAILING  (  +  Q) 
SEGMENTS  (49  -  64) 


Figure  71.-  Repeatability  of  Heat  Transfer  Data 
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Figure  72.-  Effect  of  Low  Rotation  on  Heat  Transfer 


HEAT  TRANSFER  RATIO,  Nu/Nu 


10.4  PRESSURE  LOSS  MEASUREMENTS 


Measurements  were  obtained  during  all  heat  transfer  tests  to  determine  the 
pressure  drop  through  the  serpentine  coolant  passage  (Figure  8).  The  measure¬ 
ments  for  a  nonrotating  flow  condition  at  a  lower-than-standard  flow  pressure 
showed  that  the  prinicipal  pressure  drop  occurred  in  the  180®  turns  connecting 
the  four  straight  test  section  segments.  The  uncertainty  in  the  pressure 
measurements  and  the  low  dynamic  head  of  the  flow  at  the  standard  flow  pressure 
1.01325  X  lO^N/m^  (10  atm)  precluded  the  acquisition  of  data  which  could 
yield  "benchmark  quality"  results  for  both  the  rotating  and  nonrotating  tests. 
Typical  results  from  these  data  indicate  a  variation  in  the  pressure 
distribution;  however,  the  overall  pressure  drop  from  the  inlet  to  the  exit  of 
the  model  indicated  small  effects  of  rotation.  Following  is  a  discussion  of 
the  instrumentation,  estimated  accuracy,  data  analysis  and  typical  results 
from  the  pressure  distribution  tests. 

10.4.1  Instrumentation 


The  pressure  tap  locations  are  located  in  16  sidewall  test  surface  elements  as 
shown  in  Figure  8.  The  pressure  measurements  were  obtained  with  a  Scanivalve 
Model  Z0C14  (zero,  operate,  calibrate)  differential  electronic  pressure 
scanner.  The  pressure  measurement  equipment  was  located  on  the  rotating  arm  at 
a  radius  of  approximately  30.5  cm  (12").  The  pressure  scanner  was  encased  in  a 
thermal  control  unit  (Scanivalve  Model  ZOC  TCU).  The  thermal  control  unit  is 
specified  to  maintain  a  uniform  temperature  to  ±2®C  (+3.6®F)in  a  stationary 
environment.  The  differential  pressure  transducers  were  referenced  to  pressure 
tap  #1  (Figure  8).  The  power  to  and  the  signals  from  the  electronic  pressure 
scanner  and  thermal  control  unit  were  transmitted  through  sliprings  for  the 
rotating  tests. 

10.4.2  Estimate  of  Pressure  Measurement  Accuracy 

Although  the  Z0C14  electronic  pressure  scanner  was  rated  to  be  accurate  within 
+  0.08%  full  scale  of  the  6895N/m2  (1  psi),  full  scale  value,  i.e., 

5.516N/m2  (+  0.0008  psi),  the  uncertainties  associated  with  the  transducer 
temperature  variation,  the  rotating  environment  and  the  slipring  noise  reduced 
the  estimated  accuracy  and  repeatability  to  approximately  1%  full  scale  or 
68.95  N/m^  (0.01  psi).  This  accuracy  is  compatible  with  previous  UTRC 
experience  for  pressure  measurements  in  rotating  systems.  The  value  of  68.95 
N/m2  (0.01  psi)  is  also  approximately  1.0  baseline  flow  condition. 

10.4.3  Data  Analysis 

The  pressure  drops  measured  in  the  rotating  serpentine  coolant  passage  model 
were  referenced  to  the  pressure  drop  in  the  nonrotating  model  by  (1)  subtract¬ 
ing  the  pressure  increases  due  to  work  from  model  rotation  on  the  coolant  at 
locations  radially  outward  from  pressure  tap  #1  and  (2)  adding  the  pressure 
drop  due  to  work  from  model  rotation  on  the  coolant  at  locations  radially  in¬ 
ward  from  pressure  tap  //I.  In  the  four  legged  model  (Figure  8),  this  pressure 
rise  is  of  the  order  The  pressure  increase  due  to  a  nominal  33.3®C 

(60®F)  temperature  Increase  in  the  coolant  fluid  temperature  is  approximately 
20.685  N/m^  (0.003  psi  or  0.4  Qjj^)  for  the  550  rpm  test  condition. 
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10.4.4  Results  for  0=0  and  Low  Pressure 

Pressure  drop  measurements  were  obtained  for  a  unheated,  nonrotating  test 
condition  at  the  same  Reynolds  number  as  Test  No.  1  but  at  a  lower  pressure, 
i.e.  p  =  249.599  N/m^  (36.2  psia)  compared  to  lOlTOlS/N/m^  (147.5  psia). 

This  was  done  to  increase  the  dynamic  head  of  the  flow  in  the  test  sections  to 
209.6  N/m^  (0.84  in  water  or  0.03  psi).  This  corresponds  to  about  three 
percent  of  full  scale  for  the  pressure  transducers.  The  pressure  drop  results 
are  shown  in  Figure  74  as  the  ratio  of  pressure  drop  from  the  inlet  (AP)  to 
the  dsmamic  head  of  the  air  in  the  inlet  test  section  (Qin)*  Note  the  pres¬ 
sure  drop  in  the  straight  sections  is  negligible  (as  expected)  and  is  about 
0.8  Qin  each  turn  (also  as  expected).  The  pressure  drop  measurements  will 
become  more  meaningful  in  the  rough  wall  experiments  where  the  pressure  drop 
in  the  straight  sections  will  increase  by  a  factor  of  10  and  become  measurable. 


m  =  0.0056  kg/sec  (0.0123  Ibm/sec) 

£2  =  0  rpm 

P  in  =  249599  N/m2  (36.2  psia) 

AT=0“C  (0“F) 

V  =11 .98  m/sec  (39.3  ft./sec) 

Re  =  25045 

Q  in  =  209.6  N/m^  (0.843  in  HgO  =  0.0304  psi) 

A 

4  ^  ^ 

a  A 

A 

A 

A 

A 

Jl. 

Q  I  OUTWARD  I  FIRST  I  INWARD  I  SECOND  I  OUTWARD  I  THIRD  I  INWARD 

1  STRAIGHT  I  TURN  |  STRAIGHT  |  TURN  |  STRAIGHT  |  TURN  |  STRAIGHT 

G  -  GUARD 


Figure  74.-  Pressure  Distribution  in  Nonrotating  Model  at  the  Low  Model 
Pressure  Condition. 


10.4.5  Results  for  Nonrotating  Tests 

Pressure  drops  were  measured  in  the  nonrotating  model  for  the  operating  flow 
conditions  of  Tests  1,  2  and  3  (see  Table  I).  The  pressure  drop  in  the  model 
with  respect  to  the  pressure  at  Tap  No.  1  is  shown  in  Figure  75.  The  pres¬ 
sure  drop  distributions  measured  for  the  two  higher  Reynolds  numbers  were 
similar  to  that  shown  in  Figure  74  for  a  lower  pressure  flow  condition.  Note 
that  there  is  no  measurable  pressure  drop  in  the  straight  sections  and  a 
pressure  drop  of  approximately  0.8  to  1.0  at  each  of  the  of  the  180® 
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Figure  75.-  Effect  of  Reynolds  Number  on  Pressure  Distribution  in  Nonrotating 
Model . 

turns.  As  noted  previously,  the  uncertainty  for  the  measurements  for  Test  No. 

1  are  +  1.0  Qin*  This  uncertainty  decreases  to  0.25  Qjn  Test  No.  3 
with  Re  equal  approximately  50,000  and  Increases  to  4.0  Test  No.  2 

with  Re  equal  approximately  12,500. 

10.4.6  Results  for  Rotating  Tests 

Pressure  drops  were  measured  for  rotating  Test  Nos.  4,  8  and  114  with  Re 
approximately  equal  to  25,000.  These  results  for  rotating  flow  conditions  are 
compared  with  results  for  Test  No.  1.  For  the  first  leg  of  the  coolant  passage, 
the  pressure  drops  for  all  test  conditions  are  approximately  equal  within  the 
measurement  uncertainty.  The  pressure  drop  increases  with  increasing  rotation 
rate  for  flow  in  the  second  and  third  legs.  However,  the  pressures  near  the 
coolant  passage  exit  (pressure  measurement  location  16)  are  approximately 
equal,  i.e.  ±  0.3  Qjn*  Within  the  assiamptions  made  for  the  analysis  of  the 
data,  there  is  no  explanation  for  the  large  (i.e. AP/Q^n  =  3)  variations  in 
the  pressure  distributions  shown  in  Figure  76.  However,  the  translation  of  the 
local  pressures  measured  by  the  ZOC  scanner  at  a  radius  of  approximately  30.5 
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Figure  76.-  Effect  of  Rotation  on  Pressure  Distribution  in  Model. 


cm  (12")  to  the  model  at  radii  of  68.6  cm  (27")  to  81.3  cm  (32")  requires  an 
assumption  of  the  temperature  in  the  pressure  line.  A  variation  of  temperature 
equal  to  5.6®C  (10®F)  in  the  16  lines  leading  from  the  model  to  the  ZOC 
electronic  pressure  scanner  could  account  for  the  variation  of  AP/Qjn  equal 
3.0  measured  in  the  second  leg  with  rotation. 

10.5  Model  Geometry  Information 

Table  V  gives  the  test  model  heat  transfer  geometry  information  as  a  function 
of  the  local  test  segment  heater  location  (1-64).  Table  VI  gives  the  test 
model  pressure  tap  geometry  information  as  a  function  of  pressure  tap  location 
(1-16).  The  List  of  Symbols  contains  the  colxjmn  variable  definitions. 
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TABLE  V 


MODEL  HEAT  TRANSFER  GEOMETRY  INFORMATION 


Heater  d 


A 


CA 


RADH 


S 


X 


Segment 

2 

c  2 

2 

No. 

cm 

(  in) 

cm 

(  in  ) 

cm 

1 

1.3155 

0.5179 

1.5923 

0.2468 

3.7258 

o 

C. 

1.3155 

0.5179 

1.5923 

0.2468 

5 . 2484 

3 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

4 

1.3155 

0.5179 

1.5923 

0.2465 

5.2484 

5 

1.4323 

0.5639 

1.9555 

0.3031 

2.5116 

6 

1.4323 

0.5639 

1.9555 

0.3031 

2.4277 

7 

1.4323 

0.5639 

1.9555 

0.3031 

2.4277 

8 

1.4323 

0.5639 

1 .9555 

0.3031 

2.5116 

9 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

10 

1.3155 

0.5179 

1.5923 

0 . 2468 

5 . 2484 

n 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

12 

1.5155 

0.5179 

1.5925 

0.2465 

5.2484 

13 

1.3155 

0.5179 

1.5923 

0.2468 

5 . 2484 

14 

1.3155 

0.5179 

1.5923 

0.2468 

5 . 2484 

15 

1.4323 

0.5639 

1.9555 

0.3031 

2.5116 

16 

1.4323 

0.5639 

1.9555 

0.3031 

2.4277 

17 

1.4323 

0.5639 

1.9555 

0.3031 

2.4277 

18 

1.4323 

0.5639 

1.9555 

0.3031 

2. 5116 

19 

1.3155 

0.5179 

1.5923 

0.2468 

3.7258 

20 

1.3155 

0.5179 

1.5923 

0.2468 

5 . 2484 

21 

1.3155 

0.5179 

1.5923 

0.2468 

5 . 2484 

22 

1.3155 

0.5179 

1.5923 

0.2468 

5 . 2484 

23 

1.3155 

0.5179 

1.5923 

0.2468 

5 . 2484 

24 

1 .  ail's 

0.5179 

1.5923 

0.2468 

5.2484 

25 

1.3155 

u!5179 

1.5923 

0 . 2468 

5 . 2484 

26 

1.4323 

0.5639 

1.9555 

0.3031 

2.5116 

27 

1.4323 

0.5639 

1.9555 

0.3031 

2.4277 

28 

1.4323 

0.5639 

1.9555 

0.3031 

2 . 4277 

29 

1.4323 

0.5639 

1,9555 

0.3031 

2.5116 

30 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

31 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

32 

1.3155 

0.5179 

1 . 5923 

0.2468 

5 . 2484 

33 

1.3155 

0.5179 

1.5923 

0.2468 

3.7258 

34 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

35 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

36 

1.3155 

0.5179 

1.5923 

0.2468 

5 . 2484 

37 

1.4323 

0.5639 

1.9555 

0.3031 

5.1877 

38 

1.4323 

0.5639 

1.9555 

0.3031 

5.1877 

39 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

40 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

41 

1.3155 

0.5179 

1 . 5923 

0.2468 

5 . 2484 

42 

1.4323 

0.5639 

1.9555 

0.3031 

5.1877 

43 

1.4323 

0.5639 

1.9555 

0.3031 

5.1877 

44 

1.3155 

0.5179 

1 . 5923 

0.2468 

5.2484 

45 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

46 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

47 

1.4323 

0.5639 

1.9555 

0.3031 

5.1877 

48 

1.4323 

0.5639 

1.9555 

0.3031 

5.1877 

49 

1.3155 

0.5179 

1.5923 

0.2468 

3.7258 

50 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

51 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

52 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

53 

1.4323 

0.5639 

1.9555 

0.3031 

5. 1877 

54 

1  .4323 

0.5639 

1.9555 

0.3031 

5. 1877 

55 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

56 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

57 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

56 

1.4323 

0.5639 

1.9555 

0.3051 

5.1877 

59 

1 .4323 

0.5639 

1.9555 

0.3031 

5.1877 

60 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

61 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

62 

1.3155 

0.5179 

1.5923 

0.2468 

5.2484 

63 

1.4323 

0.5639 

1.9555 

0.3031 

5.1877 

64 

1.4323 

0.5639 

1.9555 

0.3031 

5.1877 

2 
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0 . 6So 
2.4C7 

4.407 

6.407 
7.734 
8.387 
9.041 
9.694 
1  000 
3!o66 

5.000 
1 . 000 
3.000 
5.000 
6.327 
6.980 
7.634 
8.287 
0.688 

2.407 

4.407 

6.407 
1.000 
3 . 000 
5.000 
6.327 
6.980 
7.634 
8.287 
1.000 
3.000 
5.000 
0.688 

2.407 

4.407 

6.407 
8.061 
9.368 
1.000 
3.000 
5.000 

6.654 

7.961 
1.000 
3.000 
5.000 

6.654 

7.961 

0.688 

2.407 

4.407 

6.407 
8.061 
9.368 
1.000 
3.000 
:  .  OOG 

6.654 

7.961 
1 .000 
3.000 
5.000 

6.654 

7.961 
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TABLE  VI 


MODEL  PRESSURE  TAP  GEOMETRY  INFORMATION 


RADP  S 


Pressure 


Tap  No 

.  cm 

(in) 

cm 

( in) 

cm 

(in) 

1 

11.582 

4.560 

3.139 

1.236 

3.139 

1.236 

2 

16.662 

6.560 

8.219 

3.236 

8.219 

3.236 

3 

21.742 

8.560 

13.299 

5.236 

13.299 

5.236 

4 

26.822 

10.560 

18.379 

7.236 

18.379 

7.236 

5 

30.175 

11.880 

22.570 

8.886 

22.570 

8.886 

6 

26.822 

10.560 

25.491 

10.036 

0.038 

0.015 

7 

21.742 

8.560 

30.571 

12.036 

5.118 

2.015 

8 

16.662 

6.560 

35.651 

14.036 

10.198 

4.015 

9 

11.430 

4.500 

40.909 

16.106 

15.456 

6.085 

10 

8.890 

3.500 

44.770 

17.626 

19.317 

7.605 

11 

11.430 

4.500 

46.904 

18.466 

21.450 

8.445 

12 

12.217 

4.810 

47.691 

18.776 

0.358 

0.141 

13 

17.297 

6.810 

52.771 

20.776 

5.438 

2.141 

14 

22.377 

8.810 

57.851 

22.776 

10.518 

4.141 

15 

30.175 

11.880 

66.487 

26.176 

19.154 

7.541 

16 

11.582 

4.560 

84.673 

33.336 

15.674 

6.171 
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LIST  OF  SYMBOLS 


Ac 

BFMA 

Cf 

CA 

CFMA 

d 

D 

E 

FMA 

Gr 

Grx 

Gr/Re2 

Gr^/Re^ 


h 

I 

I 

j 

k 

k 

L 

L/d 


Cross  sectional  area  of  the  passage,  cm^  (in^) 

Flow  meter  (B)  reading,  kg/sec  (Ibm/sec) 

Skin  friction  coefficient  (pv2/2),  dimensionless 

Local  heated  test  surface  segment  (heater  segment)  surface 
area,  cm^  (in^) 

Flow  meter  (C)  reading,  kg/sec  (Ibm/sec) 

Local  hydraulic  diameter  of  the  passage,  cm  (in) 

Channel  depth  (Moore  1967),  cm  (in) 

Heater  element  supply  voltage,  volts 
Flow  meter  reading,  kg/sec  (Ibm/sec) 

Local  rotational  Grashof  number  based  on  hydraulic  diameter, 

_  P^,)(P^j/p2) j  dimensionless 

Local  rotational  Grashof  number  based  on  streamwise  distance  X, 
(n^  rx3)(P^,  _  Pv)(^b/p^)»  dimensionless 

Local  buoyancy  parameter  based  on  hydraulic  diameter, 

((P]3  -  P^)/P\y)(.  Cld/V) ,  dimensionless 

Lpcal  buoyancy  parameter  based  on  stre2unwise  distance  X, 

((^b  -  ^w)/  Pb)<OR/V)(nd/V)(X/d),  dimensionless 

Local  heat  transfer  coefficient,  Qnet'^(^A(T^  -  T^,), 
W/(m2c)(Btu/(hrf t2F) ) 

Unit  coordinate  direction,  (Suo  1980,  see  Section  10.1), 
dimensionless 

Heater  segment  current,  amps 

Unit  coordinate  direction,  (Suo  1980,  see  Section  10.1), 
dimensionless 

Thermal  conductivity,  W/(mC)  (Btu/(hrftF) ) 

Unit  coordinate  direction,  (Suo  1980,  see  Section  10.1), 
dimensionless 

Channel  width  (Moore  1967),  cm  (in) 

Channel  aspect  ratio  (Moore  1967),  dimensionless 
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m 

n 

Nu 

Nux 

Nuq 

Nuoo 

P 

Pin 

PMin 

Pm 

P* 

AP 

AP/Qin 

Pr 

Qin 

Qnet 

QALOSS 

r 

R 

R 

R/d 


Mass  flow  rate,  kg/sec  (Ibm/sec) 

Revolutions  per  minute  (Ito  and  Nambu  1970),  rpm 

Local  Nusselt  number  based  on  hydraulic  diameter,  hd/k, 
dimensionless 

Local  Nusselt  number  based  on  distance  from  the  channel 
entrance,  hX/k,  dimensionless 

Local  Nusselt  number  based  on  hydraulic  diameter  for  the 
stationary  (non-rotating)  test  condition,  dimensionless 

Local  Nusselt  number  for  fully  developed  turbulent  flow  in  a 
square  duct  0.0176(Re0-®) ,  dimensionless 

Pressure,  (local  pressure  (Suo  1980)  see  Section  10.1), 

N/m^  (psi) 

Model  Pressure  at  tap  #1,  N/m^  (psi) 

Inlet  model  pressure,  N/m2(psi) 

Average  channel  pressure,  (Suo  1980,  see  Section  10,1), 

N/m2  (psi) 

Non-dimensional  pressure 

Pressure  difference  between  the  local  static  pressure  tap 

readings  (corrected  for  pumping  effects  to  the  radius  of  model 

pressure  tap  #1)  and  the  model  pressure  tap  #1,  N/m^  (psi) 

Pressure  difference  between  the  local  static  pressure  tap 

readings  (corrected  for  pumping  effects  to  the  radius  of  model 

pressure  tap  #1)  and  the  model  pressure  tap  #1,  normalized  by 
the  inlet  dynamic  pressure,  dimensionless 

Prandtl  number,  dimensionless 

Channel  inlet  dynamic  pressure,  N/m^  (psi) 

Net  heat  rate  input  into  the  flow  field  from  the  heater 
segment,  W  (Btu/hr) 

Heat  flux  loss  from  the  heater  elements,  W/m^  (Btu/(hrft2) ) 

Radius  vector  (Suo  1980,  see  Section  10.1),  cm  (in) 

Radial  distance  from  the  axis  of  rotation,  RADH  +  RBASE  for 
the  heater  segments,  RADP  +  RBASE  for  the  pressure  tap 
locations,  cm  (in) 

Average  model  radius  from  the  axis  of  rotation,  cm  (in) 
Non-dimensional  radial  location 
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R/d 

Rax 

RADH 

RADP 

RBASE 

Re 

Rex 

Ro 

S 

STD 

t 

t' 

Tb 

Tw 

At 

At/t 

(  AT/T)in 

TMin 

u'i 


Non-dimensional  radial  location  based  on  the  average  model 
radius 

Rotational  Rayleigh  number  based  on  distance  from  the  channel 
entrance,  (  AT/(i;2Tj3))Pr,  dimensionless 

Radial  location  to  the  center  of  a  heater  segment  relative  to 
the  base  of  the  model  (RBASE).  cm  (in) 

Radial  location  of  wall  pressure  tap  relative  to  the  base  of 
the  model  (RBASE).  cm  (in) 

Distance  from  the  axis  of  rotation  to  the  base  of  the  model. 
The  "base  of  the  model"  was  an  arbitrarily  defined  reference 
point.  RBASE  equals  46.673  cm  (18.375_in)  for  R/d  =  49,  and 
RBASE  equals  24.130  cm  (9.50  in)  for  R/d  =  33. 

Reynolds  number  based  on  hydraulic  diameter,  md/(Ac/u  ), 
dimensionless 

Reynolds  number  based  on  distance  from  the  channel  entrance, 

P  VX/p  ,  dimensionless 

Rotation  number,  fi d/V,  dimensionless 

Distance  measured  from  the  leading  edge  of  the  first  channel 
heater  segment,  cm  (in) 

Standard  rotating  baseline  flow  conditions.  Re  =  25000, 

Ro  =  0.24,  (AP/p)in  =  (AT/T)i„  =  0.13,  R/d  =  49,  a  =  0, 
dimensionless 

Time,  (Suo  1980,  see  section  10.1),  sec 
Non-dimensional  time 

Local  bulk  temperature  of  the  fluid,  K  (R) 

Local  wall  temperature  of  the  channel  heater  segment,  K  (R) 

Temperature  difference,  (T^  -  T^j),  C  (F) 

Local  temperature  ratio,  (T^,  -  Ti3)/Tw,  dimensionless 

Temperature  ratio,  (T^  -  T|})/T^,,  where  T^,  is  evaluated 
at  the  inlet  of  the  model,  dimensionless. 

Note;  (AT/T)in  =  (AP/P)in 

Surface  temperature  at  the  model  inlet,  K  (R) 

Local  fluctuating  velocity,  m/sec  (ft/sec) 
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Ui 

u 

u/U 

V 

V 
V' 
Vin 
Vm 

VB 

VC 

VD 

X 

X’ 

X 

X/d 

y 

y 

z 
z ' 

<x 

6 

6* 

P 

P 


Time  average  local  velocity,  m/sec  (ft/sec) 

Average  channel  velocity,  m/sec  (ft/sec) 

Non-dimensional  velocity  ratio  (Moore  1967),  dimensionless 

Local  average  channel  velocity,  m/sec  (ft/sec) 

Local  velocity,  (Suo  1980,  see  Section  10.1),  m/sec  (ft/sec) 

Non-dimensional  velocity 

Inlet  channel  velocity,  m/sec  (ft/sec) 

Average  channel  velocity,  (Suo  1980,  see  Section  10.1), 
m/sec  (ft/sec) 

Thermocouple  voltage,  volts 

Thermocouple  voltage,  volts 

Thermocouple  voltage,  volts 

Coordinate  direction  (Suo  1980,  see  Section  10.1),  cm  (in) 
Non-dimensional  coordinate  direction 

Distance  measured  from  the  leading  edge  of  the  first  channel 
heater  segment  or  the  end  of  the  last  turn  region,  cm  (in) 

Non-dimensional  distance  measured  from  the  leading  edge  of  the 
first  channel  heater  segment  or  end  of  the  last  turn  region. 

Coordinate  direction  (Suo  1980,  see  Section  10.1),  cm  (in) 

Non-dimensional  coordinate  direction 

Coordinate  direction  (Suo  1980,  see  Section  10.1),  cm  (in) 
Non-dimensional  coordinate  direction 

Model  test  section  orientation  angle  with  respect  to  the  axis 
of  rotation,  deg 

Boundary  layer  thickness,  cm  (in) 

Displacement  thickness,  cm  (in) 

Absolute  viscosity,  kg/jnsec  (Ibm/ftsec) 

Kinematic  viscosity,  m^/sec  (ft^/sec) 

Local  density,  (Suo  1980,  see  Section  10.1), 
kg/m^  (Ibm/ft^) 
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Local  channel  density  based  on  local  bulk  temperature, 
kg/m^  (Ibm/ft^) 

Pjj  Average  channel  density,  (Suo  1980,  see  Section  10.1), 

kg/m^Clbm/ft^) 

Local  channel  density  based  on  local  wall  surface  temperature, 
kg/m^Clbm/ft^) 


P’ 

A  p/p 
(AP/P)in 


n 

(nR/V)(AT/T) 


Non-dimensional  density 

Local  density  ratio,  (  -  Py)/  dimensionless 

Density  ratio,  (P^  -  Py)  P^,  where  P^  is  evaluated  at 
the  inlet  of  the  model,  dimensionless. 

Note;  (AP/P)in  =  (AT/T)in 

Wall  shear  stress,  N/m^  (psi) 

Rotation  rate,  radian/sec 

Centripetal  buoyancy  parameter,  dimensionless 


Subscripts: 


b 

c 

f 

i 

in 

m 

net 

w 

X 

00 

0 


Bulk  property 

Cross  section 

Skin  friction 

At  i'th  location 

Inlet  to  model 

Average  channel  property 

Net  input  to  the  fluid 

Wall  surface  property 

Based  on  streamwise  location  X 

Fully  developed  turbulent  smooth  tube  flow 

Stationary  (non-rotating)  condition 


Superscripts: 


Vector  quantity  or  average  value 

Non-dimensional  parameter  or  fluctuating  parameter 

Dlspacement 
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